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EXECUTIVE  SUMMARY 

Mechanically  induced  residual  stresses  are  an  inherent  feature  of  permanently 
deformed  metal-matrix  composites  such  as  Al-SiC.  The  presence  of  hard  SiC 
reinforcements  promotes  non-homogeneous  flow  in  the  matrix,  resulting  in  residual 
stresses  when  the  applied  loads  are  removed.  These  residual  stresses  are  potentially 
important  because  of  the  effect  on  cavitation  and  other  damage  mechanisms  leading  to 
composite  failure.  We  have  attempted  to  simulate  the  residual  stresses  in  the  composite 
by  solving  appropriate  boundary  value  problems  using  the  finite  element  method.  Using 
this  approach,  we  have  compared  predicted  average  residual  elastic  strains  with  those 
measured  by  neutron  diffraction  for  Al-SiC  composites  subjected  to  different  loading 
histories.  The  composite  studied  was  a  2009  aluminum  alloy  reinforced  with  15  vol  % 
SiC  whiskers.  Samples  were  subjected  to  various  levels  of  plastic  strain  and  unloaded 
before  measuring  peak  shifts  in  neutron  diffraction  patterns  to  determine  residual  strains. 
Residual  strains  were  also  calculated  using  a  unit  cell  model  and  a  modified  Eshelby 
method.  The  Eshelby  method  predicted  a  linear  dependence  of  residual  elastic  strains  on 
the  imposed  deformation,  in  disagreement  with  experimental  results.  However,  the  finite 
element  calculations  produced  good  quantitative  agreement  with  experimental 
measurements  and  indicated  that  inclusion  geometry  can  significantly  affect  the  values  of 
even  average  field  quantities.  Interestingly,  the  thermally  induced  residual  strains  (from 
quenching)  were  quickly  dominated  by  the  mechanically  induced  residual  strains,  even 
for  relatively  small  amounts  of  composite  plastic  deformation.  The  excellent  agreement 
between  predicted  and  measured  residual  strains  suggest  that  it  is  possible  to  predict 
microstructural-level  stresses  with  unit  cell,  continuum-based  models  using  the  finite 
element  method.  Future  work  will  address  residual  stresses  in  particulate  composites  and 
the  effects  of  more  complex  loading  paths. 

We  have  also  been  engaged  in  studies  of  environmental  effects  on  creep  damage 
mechanisms  in  ceramic  matrix  composites.  The  benefits  of  whisker  reinforcement  of 
monolithic  ceramics  with  regard  to  fracture  toughness  is  well-established,  and  because  of 
the  inherent  thermal  and  chemical  stability,  it  is  important  to  determine  to  what  extent 
these  benefits  extend  to  high-temperature  mechanical  properties.  Samples  of  SiC 
whisker-reinforced  alumina  were  deformed  under  constant  load  at  12()0-1400°C  in 
compression.  Parallel  experiments  were  conducted  in  air  and  in  inert  atmosphere  in 
order  to  isolate  the  effects  of  test  ambient  Although  the  stress  dependence  of  the  strain 
rate  was  independent  of  test  ambient,  the  strain  rates  in  the  aerobic  ambient  were  2-5 
times  greater.  Microscopic  observations  of  crept  samples  revealed  that  glass  phase 
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accumulated  at  grain  boundary-inteiface  junctions  and  facilitated  cavitation,  despite  the 
macroscopic  state  of  uniaxial  compression.  The  glass  phase  was  generated  primarily  by 
thermal  oxidation  of  SiC  exposed  at  the  specimen  surface,  which  subsequently  penetrated 
the  composite  along  internal  boundaries.  Susceptibility  to  oxidation  is  likely  to  pose  a 
problem  for  any  carbide  reinforced  ceramic,  especially  for  tensile  loadings,  unless  the 
reaction  can  be  controlled  through  appropriate  protective  coatings. 

Presently,  we  are  investigating  the  effects  of  whisker  volume  fraction  on  creep 
response,  and  comparing  the  creep  response  in  compression  versus  flexure  using  an  inert 
test  ambient  Preliminary  results  indicate  that  the  decrement  associated  with  adding  a 
given  vol%  of  SiC  reinforcement  decreases  as  the  volume  fraction  increases,  and  that 
there  is  little  improvement  in  strain  rate  for  volume  fractions  greater  than  ~  30%. 
However,  it  is  not  clear  why  this  is  so,  or  how  the  aspect  ratio  of  the  fibers  might  affect 
the  observed  dependence.  Experiments  and  calculations  are  presently  in  progress  to 
investigate  these  issues. 

Also  included  in  this  report  is  publications  describing  the  microstructural 
evolution  and  micromechanical  modeling  of  deformation  in  2124  Al-SiC^  composites. 

This  paper  is  the  result  of  collaborative  research  involving  several  members  of  the 
engineering  faculty  at  Brown.  In  this  material,  the  yield  strength  was  relatively 
insensitive  to  aging  time  despite  a  substantial  change  in  the  yield  strength  of  the  control 
alloy.  FEM  calculations  predicted  the  qualitative  trends  observed,  although  the 
complexity  of  the  microstructurc  requires  greater  understanding  than  is  currently 
available.  Microstructural  observations  of  damage  mechanisms  occurring  during  tensile 
deformation  at  20°C  and  300°C  indicated  that  cavitation  at  fiber  ends  was  important  at 
both  temperatures.  Cavitation  was  localized  near  the  fracture  surface  at  room 
temperature,  but  occurred  through  the  gage  length  at  high  temperatures,  suggesting  that 
cavitation  may  play  an  important  role  during  creep  rupture  of  these  and  similar 
composites.  This  issue  will  be  the  subject  of  investigation  in  the  coming  year. 
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A  SiC-whisker-reinforced  alumina  composite  was  crept  in 
compression  at  1200°  to  1400°C  in  an  air  ambient  and  in 
nitrogen.  The  data  were  described  by  a  power-iaw-type  con¬ 
stitutive  relation.  The  measured  value  of  the  stress  exponent 
was  n  =  1  at  1200°C  and  n  =  3  at  1300°  and  1400°C  in  both 
ambients.  TEM  observations  were  correlated  with  the  meas¬ 
ured  creep  response  to  determine  active  deformation  mecha¬ 
nisms.  Values  of  r  =  1  were  associated  with  diffusional  creep 
and  unaccommodated  grain-boundary  sliding,  while  values 
of  R  =  3  were  associated  with  increased  microstructural 
damage  in  the  form  of  cavities.  Experiments  conducted  in 
circulated  air  resulted  in  higher  creep  rates  than  comparable 
experiments  in  nitrogen.  The  accelerated  creep  rates  were 
caused  by  the  thermal  oxidation  of  SiC  and  the  resultant 
formation  of  a  vitreous  phase  along  composite  interfaces. 
The  giassy  phase  facilitated  cavitation,  weakened  interfaces, 
and  enhanced  boundary  diffusion.  [Key  words:  creep,  alu¬ 
mina,  composites,  silicon  carbide,  oxidation.] 

I.  Introduction 

The  development  of  whisker-reinforced  ceramic  com¬ 
posites  has  created  new  applications  for  ceramics  because 
of  improvements  in  fracture  toughness  and  increased  resis¬ 
tance  to  thermal  shock,  erosion,  and  creep  when  compared  to 
monolithic  materials,'  For  example,  alumina  reinforced  with 
SiC  whiskers  is  presently  used  for  cutting  tool  inserts,  valves, 
pump  components,  and  extrusion  dies.  In  addition,  whisker- 
reinforced  composites  of  silicon  nitride,  mullite,  and  molyb¬ 
denum  disilicide  have  shown  similar  improvements  in 
mechanical  properties  and  have  potential  for  use  in  a  variety 
of  structural  applications. '■° 

Because  of  the  inherent  thermal  and  chemical  stability  of 
ceramics  and  the  benefits  of  whisker  reinforcement,  it  is  natu¬ 
ral  to  consider  the  mechanical  behavior  of  ceramic  com¬ 
posites  at  high  temperatures.  Use  of  ceramic  materials  in 
high-temperature  structural  applications  is  often  limited  by 
poor  creep  resistance.  However,  several  recent  studies  have 
shown  that  whisker  reinforcement  can  lead  to  marked  im¬ 
provements  in  creep  resistance. '■'*  Under  creep  conditions, 
the  mechanical  response  of  ceramic  composites  is  strongly  af¬ 
fected  by  the  presence  of  the  reinforcing  phase,  the  strength 
of  the  interfacial  bond  between  the  fiber  and  the  matrix,  and 
the  chemical  stability  of  the  interface.  Oxidation  reactions 
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which  occur  during  creep  can  degrade  the  interface  proper¬ 
ties,  leading  to  microstructural  damage  that  strongly  affects 
the  creep  response,"  Thus,  before  ceramic  composites  can  be 
reliably  used  in  structural  applications  at  high  temperatures,  it 
is  necessary  to  understand  the  effect  of  test  ambients  on  the 
mechanisms  of  damage  and  deformation  that  occur  during 
creep.  The  objective  of  the  present  study  is  to  determine  the 
effects  of  temperature,  stress,  and  test  ambient  on  the  com¬ 
pressive  creep  response  of  a  whisker-reinforced  ceramic  com¬ 
posite  and  to  correlate  these  findings  with  transmission 
electron  microscope  (TEM)  observations  of  crept  samples  in 
order  to  determine  mechanisms  of  deformation  and  damage. 

A  widely  studied  ceramic  composite  used  in  cutting  tool 
applications  consists  of  an  alumina  matrix  reinforced  with 
SiC  whiskers.'"’ The  composites  are  typically  fabricated 
by  hot-pressing  high-purity  alumina  powder  blended  with  SiC 
whiskers  in  volume  fractions  ranging  from  0.15  to  0.35. 
Whiskers  typically  have  aspect  ratios  of  5  to  10  and  tend  to  lie 
perpendicular  to  the  hot-pressing  axis  (HPA).  The  composite 
microstructure  is  almost  completely  free  of  residual  glass 
phase."  and  in  this  sense  constitutes  a  model  material  for 
studies  of  composite  creep  behavior.  (Ceramic  systems  such 
as  silicon  nitride  often  require  additives  for  densification.  and 
these  additives  generally  remain  in  the  composite  as  residual 
glass  phase.)  Most  studies  of  composite  creep  to  dale  have 
been  performed  in  bending  under  a  "onstant  applied  load,  and 
it  has  been  shown  that  the  creep  resistance  of  the  composites 
is  generally  far  superior  to  the  unreinforced  matrices. In 
two  of  these  studies,  creep  .-'periments  were  performed  at 
1400'  to  IbOO'C  in  air.  anr  n.icrostructural  observations  sug¬ 
gested  that  the  composi'es  deformed  by  dislocation  glide. '  ‘■ 
However,  similar  creep  tests  at  lower  temperatures  (1200'  to 
1400'C)  showed  that  the  dominant  deformation  mechanism 
was  grain-boundary  (and  interface)  sliding  resulting  from  dif¬ 
fusion.  and  that  the  sliding  was  often  unaccommodated,  caus¬ 
ing  cavitation  a'  grain-boundary-interface  junctions."  *’  All  of 
the  studies  .ii-cntioned  above  were  conducted  in  air  ambients. 
Under  th  sc  conditions,  a  thermal  oxidation  reaction  is  known 
to  occur  in  which  the  SiC  whiskers  exposed  at  the  surface  are 
oxid  /cd.'*"'  The  effect  of  test  ambient  on  ciccp  response  is 
thus  an  important  issue.  A  recent  study  showed  that  test  ambi¬ 
ents  of  argon  and  air  had  a  small  but  irreversible  effect  on  the 
compressive  creep  rates  of  similar  composites."  The  research 
described  in  the  following  sections  attemps  to  further  explore 
these  effects  and  to  identify  active  mechanisms  of  deforma¬ 
tion  and  damage. 

II.  Experimental  Procedure 

The  composite  material  selected  for  the  current  investiga¬ 
tion  was  a  SiC-whisker-reinforced  alumina.’  Microstructural 
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characterization  of  the  as-fabricated  composite  showed  an  av¬ 
erage  grain  size  of  1  to  2  pm  and  a  reinforcement  volume 
fraction  of  0.33.  Additional  microstructural  details  have  been 
described  previously.*  The  primary  processing  defects  in  the 
as-fabricated  composite  material  consisted  of  occasional 
whisker-free  channels  of  matrix.'  The  observed  structure  was 
attributed  to  a  granulation  process  often  used  to  produce  con¬ 
trolled  agglomerates  of  the  constituents  prior  to  hot-pressing.” 

Compressive  creep  experiments  were  conducted  in  an  ap¬ 
paratus  designed  to  apply  constant  stress  in  a  regulated  ambi¬ 
ent  of  continuously  circulated  dry  nitrogen  or  circulated  air.‘* 
Cylindrical  specimens  4  mm  in  diameter  and  8  mm  in  length 
were  ground  and  polished  prior  to  testing.  Creep  experiments 
were  conducted  at  1200'’  to  1400°C  under  applied  stresses  of 
25  to  250  MPa.  Samples  were  held  at  conditions  of  constant 
stress  and  temperature  until  a  constant  creep  rate  was  ob¬ 
served  for  a  least  10  h.  Plots  of  strain  rate  versus  time  were 
used  to  determine  that  a  constant  creep  rate  had  been 
reached  and  maintained.  The  stress  (or  temperature)  was  then 
increased  incrementally,  and  the  procedure  was  repieated.  An 
optical  pyrometer  sighted  on  the  silicon  carbide  platen  adja¬ 
cent  to  the  specimen  was  used  to  control  the  furnace  tem¬ 
perature.  When  the  pyrometer  registered  the  desired  test  tem¬ 
perature,  the  load  was  applied  and  the  experiment  commenced. 
The  test  temperature  was  then  maintained  to  within  TC  of 
the  preset  value.  Upon  completion  of  a  test,  samples  were 
cooled  rapidly  under  constant  stress  in  an  effort  to  retain  the 
deformed  microstructure.  The  maximum  strain  allowed  in 
each  specimen  was  39c,  thus  minimizing  nonuniform  speci¬ 
men  deformation  (“barrelling”). 

The  raicrostructures  of  as-fabricated  and  crept  specimens 
were  characterized  by  TEM.  Thin  foils  were  prepared  by  sec¬ 
tioning  a  disk  from  the  center  of  each  cylindrical  specimen. 
Disks  were  ground,  polished,  and  dimple-polished  prior  to  ion 
milling.  Thinned  specimens  were  examined  in  an  analytical 
electron  microscope  operated  at  120  kV  and  equipped  with  an 
energy  dispersive  X-ray  spectrometer  (EDS).  TEM  observa¬ 
tions  of  the  as-received  composite  revealed  a  microstructure 
that  was  virtually  free  of  residual  glass  phase.  Accumulations 
of  glass  at  triple  grain  junctions  and  at  grain-boundary-inter¬ 
face  junctions  were  rarely  observed.  Figure  1(A)  is  an  image 
of  a  typical  grain-boundary-interface  junction  in  the  as-fabri¬ 
cated  composite.  No  intragranular  pwres  or  grain-boundary 
cavities  were  observed  in  the  composite,  indicating  that  den- 
sification  was  virtually  1009(:.  Close  examination  of  the  SiC- 
AlrOj  interface  using  high-resolution  methods  revealed  a  thin 


layer  (1  nm  thick)  of  noncrystalline  material  (Fig.  1(B)).  The 
interface  film  was  presumed  to  be  SiO;,  a  native  oxide  com¬ 
monly  present  on  SiC.  However,  the  matrix  grain  boundaries 
were  ostensibly  glass-free,  within  the  resolution  of  the  dark- 
field  and  high-resolution  imaging  techniques  used  (~1  nm). 

III.  Results 

(1)  Creep  Experiments 

(A)  Inert  Ambient:  A  series  of  experiments  at  constant 
temperature  and  incrementally  increased  stress  was  pierformed 
to  determine  the  steady-state  creep  rates  of  the  composite  over 
a  wide  range  of  conditions.  The  steady-state  strain  rate  was 
ascertained  from  plots  of  strain  rate  versus  time,  as  shown 
in  Fig.  2.  After  an  initial  transient,  steady-state  creep  was 
typically  established,  and  the  associated  strain  rate  was 
determined.  Creep  rates  determined  for  a  range  of  stress' 
tempierature  conditions  are  summarized  in  Fig.  3,  a  log-log 
plot  of  the  measured  secondary  creep  rates  versus  applied 
stress  for  creep  experiments  performed  in  nitrogen  and  in  air. 
The  data  conform  to  a  power-law-type  consitutive  relation  in 
which  the  creep  rate  is  proportional  to  the  applied  stress 
raised  to  an  exponent,  n  (given  by  the  slope  of  the  lines  in 
Fig.  3).  The  value  of  the  stress  exponent  for  creep  at  1200'C 
in  nitrogen  was  approximately  1,  while  at  DOO'C.  the  stress 
exponent  was  approximately  3.5.  The  data  for  tests  at  1400''C 
appear  bilinear,  with  stress  exponents  of  n  =  1  at  low  stresses 
and  n  =  3.5  at  higher  stresses,  a  phenomenon  generally 
attributed  to  a  stress-activated  change  in  deformation  mecha¬ 
nism.  Density  measurements  were  performed  on  selected 
specimens  before  and  after  creep  experiments  to  determine 
the  cavitational  component  of  creep  strain.  At  1200'C.  cavita¬ 
tion  accounted  for  about  109?  of  the  creep  strain,  while  at 
1300°C,  the  cavitational  component  of  creep  strain  was  about 
309r.  The  total  creep  strain  in  each  specimen  was  different, 
however  (0.009  at  1200'C  and  0.015  at  1300=C). 

(B)  Oxidizing  Ambient:  A  series  of  similar  creep  experi¬ 
ments  were  conducted  in  air.  the  results  of  which  are  summa¬ 
rized  in  Fig.  3.  As  in  the  experiments  performed  in  nitrogen, 
the  data  conformed  to  a  power-law  constitutive  relation.  The 
stress  exponents  over  the  entire  range  of  applied  stresses  were 
n  =  3.3  for  13(X)°C  and  n  =  3,5  for  1400°C.  There  was  no  in¬ 
dication  of  a  bilinear  response  at  either  temperature. 

A  comparison  of  creep  data  acquired  in  nitrogen  and  in  air 
showed  that,  although  the  stress  exponents  were  similar,  the 
strain  rates  were  significantly  higher  for  tests  conducted  in  the 


Fig.  I.  Characteristics  of  as-fabricated  composite  microstructure:  (AI  a  typical  grain-boundary  interface  (GB!)  junction,  showing  an  ab¬ 
sence  of  glass  phase  (arrow);  (B)  a  high-resolution  TEM  image  of  interface  between  SiC  whisker  and  AI;Oi  matrix  show  ing  a  thin  layer  of 
amorphous  silica  (arrow). 
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Fig.  2.  Strain  rale  versus  time  curves  for  AI;Oj-SiC  composites 
subjected  to  different  stress  levels  at  1200'  and  1300'C  in  nitrogen. 


Fig.  4.  Arrhenius  plots  showing  increase  in  activation  energv  for 
creep  at  a  threshold  temperature  of  -1280  C. 


oxidizing  ambient  (Fig.  3).  The  strain  rates  increased  by  fac¬ 
tors  of  3  for  tests  performed  at  1300"C,  while  specimens  tested 
at  14(X)'C  in  air  showed  a  fivefold  increase  in  strain  rate  com¬ 
pared  to  those  crept  in  nitrogen.  At  stresses  below  50  MPa, 
the  effect  of  test  ambient  on  creep  response  diminished,  and 
at  a  stress  level  of  25  MPa  there  was  no  appreciable  effect  on 
strain  rate.  Specimens  creep-tested  in  air  at  both  temperatures 
were  discolored  (from  green  to  gray)  as  a  result  of  an  oxidation 
reaction  during  testing.  However,  specimens  tested  in  nitrogen 
showed  no  color  change. 

Activation  energies  for  secondary  creep  in  air  and  in  nitro¬ 
gen  were  obtained  from  Arrhenius  plots  (Fig.  4).  At  low  tem¬ 
perature,  creep  experiments  in  air  and  in  nitrogen  yielded 
activation  energies  of  270  and  210  kJ/mol,  respectively,  while 
at  higher  temperatures,  the  activation  energies  were  650  and 
970  kJ/mol,  respectively.  The  marked  change  in  activation 
energy  occurred  at  about  1280'C. 

It  has  been  demonstrated  that  creep  mechanisms  can  be  in¬ 
ferred  from  the  measured  value  of  the  stress  exponent  for  pure 
metals  and  monolithic  ceramic'  Stress  exponents  of  3  in 
crept  monolithic  ceramics  have  been  attributed  to  dislocation 
climb  from  Bardeen-Herring  sources,  while  stress  exponents 
of  1  have  been  attributed  to  diffusional  creep.'"  However, 
for  composites  and  most  multiphase  ceramics,  deformation 
mechanisms  cannot  be  reliably  deduced  from  creep  data  alone 
because  of  complex  interactions  between  microstructural  con¬ 
stituents  that  make  model  representations  intractable.  For 


Fig,  3.  Steadv-state  strain  rates  versus  applied  stress  for  creep  ex¬ 
periments  performed  in  nitrogen  and  in  air  at  1200  to  1400  C.  The 
data  conform  to  power-law-iype  conslituiive  relations 


this  reason,  deformation  mechanisms  in  composites  were  de¬ 
termined  by  comparing  measured  creep  response  with  direct 
observations  of  the  deformed  microstructures  using  TEM. 

(2)  TEM  Observations 

(A)  Inert  Ambient:  The  microstructures  of  samples  crept 
at  1200‘C  in  nitrogen  (which  yielded  a  stress  exponent  of  II. 
were  examined  by  TEM.  The  creep  damage  was  slight  but 
detectable.  One  of  the  features  characteristic  of  these  micro¬ 
structures  was  the  formation  of  small  (0.05  to  0.1  /um)  voids  at 
a  small  percentage  of  the  observable  grain-boundary-interface 
(GBI)  junctions  (Fig.  5(A)).  A  thin  glassy  film  was  sometimes 
observed  lining  the  interior  of  these  voids,  which  were  ob¬ 
served  only  after  creep,  and  compositional  analysis  by  EDS 
showed  that  the  phase  contained  primarily  silicon  and  alu¬ 
minum,  as  well  as  small  concentrations  of  sulfur  and  molyb¬ 
denum  (Fig.  5(B)).  A  second  microstructural  feature  that 
resulted  from  creep  at  1200'C  was  the  development  of  cor¬ 
rugated  or  faceted  alumina  grain  boundaries  (Fig  5(C)I. 
Grain-boundary  facets  in  alumina  are  often  associated  with 
boundary  migration  during  creep.*'  The  energy  of  certain 
grain  boundaries  can  be  lowered  by  the  development  of  facets, 
provided  the  resulting  facet  planes  have  sufficiently  low  en¬ 
ergy  to  offset  the  increase  in  boundary  area.  Facets  on  both 
low-angle  and  high-angle  grain  boundaries  were  observed,  and 
the  facets  were  not  detected  in  the  material  prior  to  deforma¬ 
tion.  The  facets  ranged  in  height  from  0.05  to  0.2  with  a 
mean  height  of  approximately  0.1  jzm.  while  the  step  spacing 
varied  from  0.05  to  0.8  jzm.  Gram-boundary  faceting  was  ob¬ 
served  approximately  half  as  often  as  GBI  cavities. 

The  microstruciures  of  samples  creep-tested  at  1300'  and 
14(K)'C  in  nitrogen  were  also  examined  by  TEM.  The  samples 
were  characterized  by  voids  at  GBI  junctions,  at  whisker- 
whisker  contact  areas,  and  at  alumina  triple  grain  junctions, 
as  well  as  by  corrugated  grain  boundaries.  Figure  6  shows  a 
typical  cavity  approximately  0.1  /zni  in  diameter  at  a  triple 
grain  junction.  Cavities  at  these  sites  were  typically  triangu¬ 
lar.  with  small  amounts  of  glassy  phase  at  the  corners.  The 
voids  at  GBI  junctions  were  larger  (0.05  to  0.3  jtml  and  more 
frequent  (5*^  to  lOT  of  the  junctions)  than  in  specimens  tested 
at  1200'C.  Also,  voids  at  GBI  junctions  were  approximately 
twice  as  frequent  as  voids  at  triple  grain  junctions  Voids  in 
these  specimens  often  had  a  thin  lining  of  glassy  phase,  as 
shown  in  Fig.  6,  but  approximately  one-third  of  the  voids  ap¬ 
peared  to  be  glass-free.  Distinguishing  between  glass-lined 
and  glass-free  voids  is  not  always  possible,  because  glass  is 
often  preferentially  etched  by  ion  milling  and  it  can  redis¬ 
tribute  during  cooling.  Also,  glass  films  can  be  extremely  thin 
and  thus  difficult  to  resolse,  even  by  TEM.  However,  the 
observation  of  more  widespread  cavitation  indicated  that  the 
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Fig.  S.  Structural  and  chemical  features  of  the  boundary  regions  of  the  composite  after  creep  at  1200°C  in  nitrogen.  (A)  A  cavity  at  a 
grain-boundarv-interface  (GBI)  junction.  A  thin  layer  of  glass  phase  lines  the  cavity.  (B)  An  EDS  spectrum  showing  the  composition  of  a 
typical  glass  phase  accumulation  The  elements  present  are  primarily  silicon  and  aluminum  (and  presumably  oxygen)  with  traces  of  sulfur 
and  molybdenum.  (C|  A  corrugated  alumina  grain  boundary. 


Fig.  «.  A  thin  film  of  glass  lining  a  cavity  at  an  alumina  triple 
grain  junction  after  creep  in  nitrogen  at  IBOO'C.  The  glass  phase  is 
most  visible  in  the  corners 


cavitational  component  of  strain  increased  at  the  higher  test 
temperature  and  contributed  to  the  change  in  stress  exponent 
and  activation  energy. 

Corrugated  grain  boundaries  were  also  observed  in  speci¬ 
mens  tested  at  1300"^  and  1400'C.  The  step  height  and  step 
separation  of  these  grain  boundaries  were  unchanged  w  ith  re¬ 
spect  to  specimens  tested  at  12(X)'C.  Samples  tested  in  nitro¬ 
gen  at  1400'C  also  exhibited  voids  along  grain  boundaries,  and 
in  rare  cases  the  voids  were  accompanied  by  dislocations. 
Figure  7  shows  an  example  of  typical  cavities  along  an  alu¬ 
mina  grain  boundary.  The  faceted  cavities  are  approximately 
0.2  to  0.5  Aim  in  diameter  and  apparently  free  of  glass  phase. 

(B)  Oxidizing  Ambient:  The  primary  microstructural 
difference  between  specimens  crept  in  air  and  those  tested  in 
nitrogen  was  the  abundance  of  glassy  phase  in  the  former.  The 
development  of  an  amorphous  phase  in  these  specimens  led 
to  a  slight  increase  in  the  size  and  number  of  cavities  at  triple 
grain  junctions  and  at  sites  of  whisker-whisker  contact.  Speci¬ 
mens  after  creep  exhibited  glass-filled  pockets  and  glass-lined 
cavities  (0.5  to  1.5  jim)  that  occurred  at  GBI  junctions,  triple 
grain  junctions,  and  at  regions  of  whisker-whisker  contact. 
Figure  8(A)  shows  a  cavity  at  a  GBI  junction  that  was  typical 
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Fig.  7.  Faceted  voids  along  alumina  grain  boundaries  in  specimens 
crept  at  )400'C  in  nitrogen  Voids  appeared  to  be  associated  with 
gram-boundary  facets. 
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of  specimen.s  tested  in  air  at  1300'  and  1400'C.  Note  that  the 
film  lining  the  cavity  is  thicker  (about  10  nm)  than  the  film 
shown  in  Fig.  5(A)  w'here  it  is  barely  detectable.  Nearly  all  of 
the  cavities  observed  in  these  samples  were  glass-lined  and 
assumed  shapes  that  ranged  from  triangular  (for  smaller  cavi¬ 
ties)  to  polygonal  (for  larger  ones).  The  EDS  spectrum  from 
the  secondary  phase  showed  primarily  silicon  and  aluminum, 
with  small  amounts  of  calcium  and  yttrium  (Fig  8(B)).  It  was 
presumed  to  contain  oxygen  as  well,  although  the  spectrom¬ 
eter  used  could  not  detect  oxygen.  Figure  8(C)  shows  another 
cavity  that  developed  from  an  accumulation  of  glass  pha.se  at 
a  GBl  junction.  Ligaments  extend  across  the  void,  indicating 
that  the  aluminosilicate  phase  behaved  as  a  viscous  liquid  dur¬ 
ing  creep  deformation.  Corrugated  grain  boundaries  were  also 
observed,  although  associated  dislocations  were  not,  and  the 
overall  dislocation  density  was  negligible. 

Elemental  X-ray  analysis  was  performed  on  over  20  glass 
pockets  in  these  specimens,  and  the  results  showed  distinct 
differences  in  composition.  The  two  typical  EDS  spectra  en¬ 
countered  are  shown  in  Figs.  5(B)  and  8(B).  The  spectrum  in 
Fig.  5(B)  shows  predominantly  silicon  and  aluminum,  with 
small  amounts  of  sulfur  and  molybdenum  The  second  spec¬ 
trum,  Fig.  8(B)  shows  similar  concentrations  of  silicon  and 


Fig.  8.  Microsiructural  damage  altribuied  lo  creep  al  l.W  and  MOOT  in  air  lA)  glass-lined  casiu  al  GBl  junciion  in  sample  afier  creep 
in  air  al  1.400  C.  iB'  F.DS  specirum  from  glass  phase,  show  mg  impurilies  of  Ca  and  Y.  (C)  cas  iialion  wuhin  glass  phase  accumulaied  ai  GBl 
junclion  in  specimen  creep-iesied  in  air  al  14CK'  C  Glassy  ligaments  bridge  the  separated  interface 
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aluminum,  with  small  amounts  of  calcium  and  yttrium.  This 
spectrum  was  typical  of  glass  phase  present  in  composites 
after  creep  in  air.  An  additional  constituent,  although  not  de¬ 
tectable  by  the  spectrometer,  was  oxygen,  and  the  phase  can 
be  described  as  an  aluminosilicate  glass.  High  aluminum  con¬ 
tents  were  consistently  observed,  and  the  compositions  were 
well  above  the  equilibrium  solubility  of  alumina  in  silica.  The 
compositions  (including  the  impurity  content)  were  consistent 
and  included  impurities  of  either  molybdenum  and  sulfur  in 
some  cases,  or  calcium  and  yttrium  in  others.  Glass  phase  of 
both  compositions  was  often  present  in  the  same  sample. 

TEM  foils  sectioned  near  the  surface  of  samples  crept  in 
air  revealed  that  occasionally  SiC  whiskers  below  the  surface 
were  also  oxidized.  Figure  9  shows  an  example  of  a  partially 
oxidized  whisker  located  at  a  depth  of  approximately  1  mm 
below  the  surface  of  the  creep  sample.  The  whisker  flank  is 
indented  and  a  thick  film  of  amorphous  silica  has  formed 
where  the  whisker  has  been  oxidized.  Samples  were  sectioned 
from  different  depths,  and  this  feature  generally  was  not  ob¬ 
served  near  the  center  of  the  specimen. 

IV.  Discussion 

(I)  Creep  Deformation  Mechanisms 

As  described  above,  creep  experiments  at  1200‘C  in  nitro¬ 
gen  were  characterized  by  a  stress  exponent  value  of  n  =  1, 
and  only  subtle  changes  in  the  microstructure  were  observed 
by  TEM.  The  microstructural  changes  included  small  cavities 
at  GBI  junctions  and  grain-boundary  facets.  There  was  an 
insufficient  number  of  dislocations  in  the  matrix  to  support  a 
dislocation-controlled  creep  mechanism.  Previous  creep  stud¬ 
ies  of  unreinforced  alumina  have  shown  that,  for  test  condi¬ 
tions  comparable  to  those  used  here,  the  material  deformed  by 
a  grain-boundary  diffusional  creep  mechanism  and  exhibited 
a  stress  exponent  value  of  n  =  1  in  bending  "  Furthermore. 
TEM  observations  revealed  that  crept  alumina  was  character¬ 
ized  by  unaccommodated  grain-boundary  sliding,  which  re¬ 
sulted  in  cavities  at  triple  grain  junctions,  corrugated  alumina 
grain  boundaries,  and  a  low  dislocation  density."  The  simi¬ 
larities  between  observations  reported  here  for  composite 
samples  creep-tested  in  compression  and  those  reported  pre¬ 
viously  for  monolithic  aluminas  suggested  that  the  dominant 
creep-deformation  mechanism  for  the  composite  tested  under 
these  conditions  was  grain-boundary  diffusional  creep  (Coble 
creep),  with  occasional  unaccommodated  sliding  at  interfaces 
and  grain  boundaries. 

As  shown  in  Figs.  3  and  4.  as  the  test  temperature  was  in¬ 
creased  from  1200'  to  1300‘C  in  the  nitrogen  ambient,  the 
strain  rates  showed  an  increased  dependence  on  the  applied 


Fig.  1  Cavitation  at  a  subsurface  SiC  whisker  after  creep  in  air 
Thermal  oxidation  caused  (he  indentation  in  the  whisker  flank  and 
generated  glass  phase. 


stress  and  an  increase  in  activation  energy.  An  increase  in 
stress  exponent  and/or  activation  energy  is  generally  attributed 
to  the  activation  of  a  new  deformation  mechanism(s).  For  ex¬ 
ample.  tensile  creep  tests  performed  on  siliconized  silicon 
carbide  composites  showed  that  an  increase  in  stress  exponent 
was  caused  by  a  stress-activated  increase  in  microstructural 
creep  damage.^'  At  low  stresses,  creep  deformation  was  con¬ 
trolled  by  dislocation  glide.  However,  above  a  critical  stress 
level,  cavitation  occurred  at  the  matrix-reinforcement  inter¬ 
face,  and  the  stress  exponent  increased. The  TEM  observa¬ 
tions  reported  here  show  that,  for  specimens  tested  at  1300'  to 
1400'C,  cavities  were  larger  and  more  frequent  than  at  1200'C 
and  appeared  at  additional  sites  in  the  microstructure,  includ¬ 
ing  triple  grain  junctions,  grain  boundaries,  and  regions  of 
whisker-whisker  contact  (Figs.  5-7).  The  increase  in  creep 
damage  was  associated  with  the  observed  increases  in  stress 
exponent  and  activation  energy.  There  was  insufficient  evi¬ 
dence  to  support  a  dislocation  glide  mechanism  at  either 
temperature. 

Creep  deformation  mechanisms  in  multiphase  ceramics 
cannot  always  be  reliably  inferred  from  measured  values  of 
activation  energies  or  stress  exponents.  For  example,  the  acti¬ 
vation  energy  determined  for  creep  below  1280'C  was 
210  kJ/mol  in  nitrogen  and  269  kJ/mo!  in  air.  while  at  higher 
temperatures,  the  activation  energy  was  966  kJ/mol  in  nitro¬ 
gen  and  655  kJ/mol  in  air.  These  values  are  not  entirely 
consistent  with  reported  activation  energies  for  boundary 
diffusion  in  alumina  (419  kJ/mol  for  boundary  diffusion  of 
aluminum  cations  and  380  kJ/mol  for  oxygen  anions)."’  How¬ 
ever,  the  TEM  observations  did  not  support  a  dislocation 
glide  mechanism  in  the  matrix,  leaving  diffusional  creep  as 
the  most  probable  deformation  mechanism.  Note  that  when 
both  lattice  diffusion  and  grain-boundary  diffusion  are  impor¬ 
tant.  the  activation  energy  may  not  be  unique.  Because  of  the 
extremely  small  grain  size  (1  to  2  ^m).  the  relative  contribu¬ 
tion  of  boundary  diffusion  is  expected  to  be  greater  than  that 
of  lattice  diffusion. 

The  value  of  the  stress  exponent  in  =  3)  can  be  interpreted 
in  different  ways.  Warshaw  and  Norton"’  showed  that,  for  un- 
reinforced  aluminas,  the  ceramic  creep  rates  decreased  as 
grain  size  increased,  and  the  stress  exponent  of  large-grained 
alumina  was  approximately  4.  compared  to  a  value  of  n  =  I 
for  fine-grained  material.  Jakus  reported  creep  rates  of  fine¬ 
grained  whisker-reinforced  alumina  that  were  similar  to 
creep  rates  of  large-grained  alumina.  He  reasoned  that 
whisker  reinforcements  could  bond  several  grains  into  aggre¬ 
gates  that  would  effectively  behave  as  larger  alumina  grains, 
thus  resulting  in  lower  creep  rates."  Chokshi  and  Porter,  how¬ 
ever.  concluded  that  most  of  the  applied  load  was  carried  by 
the  reinforcements  and  compared  measured  stress  exponents 
(n  =  5.2)  to  those  of  monolithic  SiC  (n  =  5.7).  Cannon  and 
Sherby  suggested  that  an  exponent  ofn  =  3  could  result  from 
a  combination  of  diffusional  creep  and  dislocation  glide  proc¬ 
esses,  for  which  the  stress  exponents  are  expected  to  be  n  =  1 
and  ft  =  5  respectively.’*’  A  similar  interpretation  might  be 
applied  to  the  data  presented  here  by  arguing  that  diffusional 
creep  of  the  matrix  should  lead  to  a  low  stress  exponent 
(n  =  1),  while  creep  of  the  SiC  reinforcements  should  lead  to 
a  higher  stress  exponent  (reportedly  n  =  5.7  at  1500'C’  )  The 
strain  rate  of  the  composite  might  then  be  expected  to  show 
an  intermediate  dependence  on  the  applied  stress  between 
that  of  unreinforced  alumina  (n  =  1)  and  that  of  pure  SiC 
(n  =  5.7).  However,  no  evidence  of  creep  deformation  was 
observed  in  the  whiskers,  nor  was  there  any  evidence  of  dis¬ 
location  glide  in  the  matrix.  These  observations  support  the 
previous  assertion  that  the  composite  deformed  by  diffusional 
creep  with  unaccommodated  sliding  (caviiation)  at  1200  C. 
while  at  higher  temperatures,  where  the  stress  exponent  and 
activation  energy  were  higher,  creep  was  characterized  by  in¬ 
creased  microstructural  damage  in  the  form  of  cavitation. 

The  data  from  creep  experiments  at  1400'C  showed  an  ap¬ 
parent  stress-activated  change  in  stress  exponent  from  n  =  1 


1246 


Journal  of  the  American  Ceramic  Society  —  Lipetzky  et  al. 


Vol.  74,  No.  6 


to  n  =  3  (Fig.  3).  Crept  microstructures  (for  which  n  =  3) 
showed  that,  in  addition  to  porosity  at  GBI  and  triple  grain 
junctions,  cavities  were  also  observed  at  grain-boundary  facets 
(Fig.  7).  When  irregularly  shaped  grain  boundaries  slide  dur¬ 
ing  creep,  tensile  stress  can  develop  perpendicular  to  the  di¬ 
rection  of  sliding,^*  causing  voids  to  nucleate.  A  plausible 
sequence  leading  to  the  microstructure  shown  in  Fig.  7  would 
thus  begin  with  the  development  of  corrugated  boundaries 
followed  by  void  nucleation.  Some  of  the  cavitated  grain 
boundaries  were  also  accompanied  by  dislocations.  However, 
the  overall  dislocation  density  was  negligible,  and  it  is  unlikely 
that  dislocation  glide  played  a  significant  role  in  creep  defor¬ 
mation  at  these  temperatures.  The  presence  of  whisker  rein¬ 
forcements  undoubtedly  inhibited  grain-boundary  sliding  by 
effectively  pinning  the  matrix  grains.  However,  at  lower  tem¬ 
peratures  and  stress  levels,  most  of  the  sliding  was  accommo¬ 
dated  by  diffusion,  while  at  higher  temperatures  and  stress 
levels  the  sliding  was  increasingly  unaccommodated  Density 
measurements  supported  this  assertion,  indicating  that  the 
cavitational  component  of  creep  strain  increased  from  about 
lO^c  of  the  total  strain  at  1200°C  to  about  3Q9c  of  the  total 
strain  at  1300'C.  It  was  concluded  that  the  observed  threefold 
increase  in  stress  exponent  of  the  composite  was  associated 
with  an  increase  in  the  cavitational  component  of  strain,  al¬ 
though  other  causes  may  have  been  involved  as  well.  Finally, 
corrugated  grain  boundaries  were  also  observed,  indicating 
that  although  the  stress  exponent  and  strain  rates  increased, 
diffusional  creep  remained  active  at  these  temperatures. 

(2)  Oxidation  Effects 

The  effects  of  the  oxidizing  ambient  on  the  creep  behavior 
of  SiC-whisker-reinforced  alumina  can  be  seen  in  Figs.  3  and 
4  and  in  the  TEM  observations  of  the  crept  microsiructures 
(Figs,  8  and  9).  Similar  stress  exponents  and  microstructures 
were  observed  for  specimens  tested  in  air  and  in  nitrogen.  For 
example,  cavities  were  detected  at  the  same  types  of  locations 
in  both  air-tested  and  nitrogen-tested  microstructures,  and 
corrugated  grain  boundaries  were  observed  in  both.  However, 
the  glassy  phase  was  far  more  prevalent  in  air-tested  samples, 
and  cavities  were  larger.  In  addition,  the  strain  rates  were 
consistently  higher  than  for  comparable  stress-temperature 
conditions  in  nitrogen.  This  behavior  was  attributed  to  the 
increased  amount  of  glassy  phase  that  developed  during  creep 
in  air.  The  presence  of  viscous  amorphous  films  along  inter¬ 
faces  undoubtedly  facilitated  grain-boundary  and  interface 
sliding,  decreased  the  interface  bond  strength,  and  decreased 
the  creep  resistance  of  the  composite.  A  layer  of  silica  glass  at 
grain  boundaries  and  interfaces  would  also  facilitate  alumina 
dissolution-reprecipitation  reactions. The  aluminosilicate 
glass  film  along  boundaries  would  permit  alumina  to  dissolve 
and  reprecipitate,  thus  allowing  grains  to  change  shape  and 
accommodate  the  applied  stress.*” 

Perhaps  the  most  important  effect  of  the  oxidizing  ambient 
on  creep  was  the  role  of  the  vitreous  phase  in  facilitating  dif¬ 
fusion  and  cavitation.  This  finding  is  consistent  with  a  previ¬ 
ous  investigation  of  creep  of  vitreous-bonded  alumina,  where 
it  was  found  that  cavities  nucleated  in  the  amorphous  phase 
that  existed  at  triple  grain  junctions,'*'  In  the  composite,  the 
formation  of  cav  ities  at  GBI  junctions  and  at  triple  grain  junc¬ 
tions  allowed  matrix  grains  to  partly  overcome  the  pinning 
effect  of  the  whiskers  and  thus  slide  without  accommodation. 
The  widespread  cavitation  enhanced  by  the  abundance  of 
glass  phase  accounted  for  the  increase  in  strain  rates  in  air. 
Thus,  it  was  concluded  that  the  air-tested  samples  deformed 
by  diffusional  creep  and  that  the  boundary  sliding  was  fre¬ 
quently  unaccommodated,  resulting  in  extensive  cavitation. 

The  source  of  the  abundant  glassy  phase  that  developed 
during  creep  in  air  was  related  to  a  thermal  oxidation  reaction 
that  occurred  during  testing.  Thermal  oxidation  experiments 
on  similar  composites  have  shown  that  the  material  undergoes 
an  oxidation  reaction  when  exposed  to  air  at  temperatures 


similar  to  those  used  in  the  present  research.'*  '*  Whiskers  ex- 
piosed  at  the  surface  oxidize,  forming  silica  glass  and  carbon, 
and,  if  temperature  (and  time)  are  sufficient,  the  silica  will 
react  with  alumina  to  form  mullite.  The  glassy  phase  that  de¬ 
veloped  in  the  air-tested  specimens  originated  primarily  from 
the  oxidation  of  SiC  whiskers.  Under  the  test  conditions,  a 
silicon-rich  amorphous  phase  formed  on  the  surface  of  the 
composite  and  subsequently  penetrated  the  composite  along 
grain  boundaries  and  interfaces,  thereby  weakening  the  ma¬ 
terial.  Furthermore,  oxygen  diffusion  rates  were  undoubtedly 
higher  along  interfaces  and  boundaries  because  of  the  pres¬ 
ence  of  glass  phase,  causing  oxidation  of  internal  (subsurface) 
whiskers  in  the  composite  as  well  (Fig.  9).  Other  possible 
sources  of  the  amorphous  phase  observed  in  crept  samples 
include  redistribution  of  residual  sintering  aids  (such  as  MgO, 
Y;03,  and  CaO)  and  native  oxide  surface  films  on  the  SiC 
whiskers.  The  amount  of  glass  arising  from  these  sources  was 
apparently  much  smaller  than  the  amount  formed  dy  namically 
by  thermal  oxidation  of  SiC  whiskers. 

The  compositions  of  the  glass  phases  shown  in  Figs  5  and  8 
were  different  from  the  composition  expected  from  SiC  oxi¬ 
dation  alone.  The  presence  of  aluminum  (Figs.  5(B)  and  8iBii 
implies  substantial  solubility  of  alumina  in  the  vitreous  silica. 
It  has  been  reported  that  a  binary  glass  can  be  formed  con¬ 
taining  up  to  10  mol'r  (16  wt9r)  alumina  and  90  moUr  silica.  ' 
However,  the  intensity  of  the  aluminum  peak  in  Figs.  5iBi 
and  8(B)  corresponds  to  far  more  than  lO'T  alumina,  the  solu¬ 
bility  limit.  While  the  characteristic  X-rays  of  silicon  will  fluo¬ 
resce  and  generate  additional  X-rays  of  aluminum  from  the 
surrounding  matrix,  it  is  unlikely  that  this  could  cause  the 
observed  Al  peak  intensity.  Furthermore,  high-resolution  im¬ 
ages  did  not  reveal  the  formation  of  microcrystallites  of  alu¬ 
mina.  silica,  or  mullite.  indicating  that  the  alumina  was 
dissolved  in  a  homogeneous  glass  phase.  However,  low-level 
impurities  such  as  Mo  and  S  (Fig.  5(B(.  or  Ca  and  Y  (Fig.  8(Bu 
were  consistently  observed  in  the  glassy  accumulations,  and  it 
has  been  shown  that  these  and  other  impurities  can  greatly 
increase  the  solubility  of  alumina  in  aluminosilicate  glass.*' 
For  example,  the  addition  of  small  amounts  of  certain  oxides, 
such  as  CaO,  enables  melts  composed  of  up  to  50  mol^r  alu¬ 
mina  in  silica  to  easily  be  cooled  to  an  amorphous  state/'  The 
presence  of  impurity  elcment(s)  such  as  those  detected  in  the 
different  glass  phase  accumulations  is  likely  to  inhibit  phase 
separation  and  emulsion  during  cooling,  resulting  in  the  ob¬ 
served  homogeneous  alumina-rich  silicate  glass. 

V.  Summary  and  Conclusions 

The  present  investigation  has  shown  that  atmospheric  ef¬ 
fects  play  a  strong  role  in  the  compressive  creep  behavior  of 
AI;0)-SiC.  composites.  At  12(X)  C.  creep  in  nitrogen  oc¬ 
curred  primarily  by  diffusion,  with  some  unaccommodated 
sliding  at  interfaces  and  grain  boundaries  At  higher  tempera¬ 
tures  and  stresses  in  nitrogen,  the  sliding  was  increasingly  un¬ 
accommodated  and  the  cavitational  component  of  creep  strain 
increased,  resulting  in  a  stronger  dependence  on  the  applied 
stress  and  a  stress  exponent  value  of  3.5.  The  activation  en¬ 
ergy  for  creep  increased  from  210  to  966  kj/mol  when  the  test 
temperature  was  increased  above  1280  C.  Diffusional  creep 
continued  to  be  an  active  process  under  all  stress  temperature 
conditions  in  nitrogen  Creep  experir^ents  conducted  in  air 
showed  essentially  the  same  stress  dependence,  although 
creep  rates  were  consistently  higher.  Oxidation  of  SiC  whis¬ 
kers  exposed  at  the  specimen  surface  led  to  the  formation 
of  an  aluminosilicate  glassy  phase  that  subsequently  pene¬ 
trated  along  grain  boundaries  and  interfaces  and  accumulated 
at  GBI  and  triple  grain  junctions  throughout  the  material 
The  formation  of  glass  phase  from  thermal  oxidation  of  SiC 
strongly  affected  the  creep  behavior.  Glass  films  weakened 
internal  boundaries  in  the  material  and  increased  grain¬ 
boundary  diffusion  rates,  which  contributed  to  accelerated 
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creep  rales.  The  contribution  of  creep  damage  to  the  defor¬ 
mation  process  was  larger  for  tests  conducted  in  air  than  for 
similar  tests  in  nitrogen. 

Finally,  it  appears  that  all  ceramic  composites  reinforced 
with  carbide  whiskers  will  suffer  similar  problems  of  thermal 
oxidation  at  high  temperatures.  Significant  improvements  in 
creep  resistance  would  be  realized  by  controlling  the  proper¬ 
ties  of  the  glass  phase  through  additives  that  would  reduce 
the  viscosity  of  the  glass  or  enhance  devitrification,  or  by  the 
use  of  protective  coatings.  Alternatively,  the  development  of 
oxidation-resistant  reinforcements  such  as  oxides,  nitrides,  or 
silicides  would  be  beneficial  to  the  advancement  of  ceramic 
matrix  composites  for  high-temperature  structural  appli¬ 
cations  Successful  solutions  would  depend  on  the  stress- 
temperature  conditions  involved  in  the  application  and  on  the 
particular  service  environment. 
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Introduction 


Mechanically  induced  residual  stresses  are  an  inherent  feature  of  permanently  deformed  metal- 
matrix  composites.  The  presence  of  hard  reinforcement  particles  promotes  non-homogeneous  flow 
in  the  material,  resulting  in  residual  str-esses  w'hen  the  applied  loads  are  removed.  Several  distinct 
approaches  have  been  used  in  attempts  to  quantify  thermally  and  mechanically  induced  residual 
stresses  in  metal-matrix  composites.  One  approach  has  been  to  assume  that  the  composite  consists  of 
a  periodic  array  of  whiskers  embedded  throughout  the  matrix  [l].  With  this  assumption,  a  boundarj' 
value  problem  can  be  formulated  on  a  single  unit  cell  and  solved  by  use  of  the  finite  element  method. 
A  second  approach  uses  modified  techniques  originally  developed  by  Eshelby  [2]  to  estimate  residual 
stresses  for  composites  reinforced  by  ellipsoidal  inclusions  [3]. 

The  main  objective  of  the  present  investigation  is  to  compare  predicted  average  residual  elastic 
strains  with  those  measured  by  neutron  diffraction  for  Al/SiC  composites  subjected  to  a  rapid  quench 
followed  by  different  loading  histories.  Similar  comparisons  of  neutron  diffraction  measurements  of 
thermally  induced  residual  strains  and  elastic  calculations,  sometimes  including  a  simple  model  for 
creep  relaxation,  were  performed  on  a  variety  of  composite  materials  by  Majumdar  et  al.  [4],  Other 
objectives  are  to  compare  the  results  of  finite  element  and  Eshelby  model  predictions  of  residual 
strains  and  to  ascerteun  the  effect  of  inclusion  shape  on  predicted  values  of  average  residual  elastic 
strains. 


Experimental 

The  composite  studied  w'as  a  2009  aluminum  alloy  ( Al-3.8%Cu-1.3%Mg-0.02%Fe)  reinforced 
with  15  vol9c  silicon  carbide  whiskers,  obtained  from  the  Advanced  Composite  Materials  Corp.  in 
Greer,  S.C.  The  2009  matrix  alloy  was  specifically  designed  for  use  in  metal-matrix  composites.  The 
samples  were  solutionized  for  1  hour  at  a  temperature  of  493° C,  quenched  in  water,  and  then  deformed 
in  uniform  tension  or  compression  to  varying  levels  of  pla.stic  strain.  In  one  case,  a  specimen  was 
strained  until  fracture.  The  specimens  were  subsequently  stored  in  liquid  nitrogen  to  prevent  recovery 
in  the  aluminum  matrix.  Approximately  8  hours  prior  to  a  neutron  diffraction  experiment,  samples 
were  returned  to  room  temperature.  To  ascertain  the  mechanical  properties  of  2009  aluminum  at 
these  temperatures,  tensile  tests  w’ere  also  performed  on  the  unreinforced  matrix  in  liquid  nitrogen 
(_196°C),  in  a  solution  of  dry  ice  and  methyl  alcohol  (  — 75°C),  and  at  room  temperature.  The 
unreinforced  alloy  received  the  identical  solutionizing  treatment  given  the  composite,  although  this 
does  not  necessarily  lead  to  the  same  microstructural  state  as  is  in  the  composite  matrix  [5]. 

The  neutron  diffraction  measurements  were  performed  on  the  Neutron  Powder  Diffractometer  at 
the  Manuel  Lujan,  Jr.  Los  Alamos  Neutron  Scattering  Center  (L.4NSCE),  located  at  the  Los  Alamos 
National  Laboratory.  This  facility  utilizes  a  high-intensity  pulsed  spallation  source.  In  addition  to  the 


coipposite  samples,  experiments  were  also  performed  on  the  unreinforced  cdloy  and  on  unrestrained 
silicon  ceirbide  whiskers  to  determine  stress-free  lattice  constants.  Specimens  were  placed  in  the 
diffractometer  with  the  deformation  axis  in  the  diffraction  plane  ajid  oriented  at  45  degrees  to  the 
incident  neutron  beam.  The  samples  were  6.35  mm  in  diameter  and  38.1  mm  in  length.  Time-of- 
flight  diffraction  patterns  were  collected  in  each  of  two  opposed  detectors  placed  at  ±90  degrees  with 
respect  to  the  incident  beam  over  an  eight  hour  (typic^)  period.  The  diffraction  vectors  for  these 
patterns  were,  respectively,  parallel  and  perpendicular  to  the  deformation  axis. 

Bragg  reflections  in  each  spectrum  were  fitted  individually  to  obtain  interplanar  spacings.  In 
addition,  lattice  parameters  were  obtained  for  the  aluminum  matrix  by  use  of  Rietveld  profile  re¬ 
finement  [6].  One  advantage  to  profile  refinement  is  that  all  the  lattice  reflections  contribute  to 
the  determination  of  the  lattice  parameter.  Once  the  interplanar  spacings  axe  obtained  from  the 
diffraction  patterns,  the  lattice  strains  are  then  given  by  the  relations 


^hkl  = 


dkkl  - 


hkl 


(1) 


where,  for  example,  dhkl  represents  matrix  lattice  spacings  in  the  composite  and  represents 
corresponding  lattice  spacings  in  the  unreinforced  alloy,  and  e/tki  is  the  calculated  matrix  strain  for 
a  given  lattice  plane.  For  the  profile  refinements,  the  strain  is  calculated  by  replacing  values  of 
interplanar  spacings  in  (1)  by  the  values  of  the  lattice  parameters. 


Theory 


Withers  et  al.  [3]  developed  an  analytical  method  that  uses  solutions  for  ellipsoidal  inclusions 
provided  by  Eshelby  [2]  to  determine  the  residual  stress  states  in  thermally  and  mechanically  deformed 
composites.  In  reference  [3]  relations  are  developed  that  express  the  inclusion  transformation  strain 
in  terms  of  the  plastic  deformation  undergone  by  the  composite.  From  the  calculated  transformation 
strain,  the  average  residual  stresses  and  strains  are  predicted.  The  chief  advantage  to  this  approach 
is  that  no  extensive  computing  is  required. 

Although  axisymmetric  unit  cell  based  finite  element  ceilculations  require  larger  amounts  of  com¬ 
puter  time  (but  still  considerably  less  than  a  full  three  dimensional  simulation),  they  can  incorporate 
more  gener^  inclusion  geometries  and  material  properties.  For  the  computations  here,  the  method¬ 
ology  of  reference  [7]  is  used.  Briefly,  the  material  is  modelled  as  a  periodic  eirray  of  hexagonal 
cylinders,  each  consisting  of  the  aluminum  alloy  with  a  silicon  carbide  fiber  embedded  in  the  cen¬ 
ter.  The  hexagonal  cylinder  can  then  be  approximated  as  a  circular  cylinder,  so  that  the  problem  is 
rendered  axisymmetric.  In  addition,  because  of  symmetry,  only  one  quadrant  of  the  cylinder  needs 
to  be  analyzed.  Both  right  circular  cylindrical  and  ellipsoidal  reinforcements  are  considered;  for  the 
calculations  here  the  aspect  ratio  of  the  unit  cell  and  of  the  reinforcement  is  tahen  to  be  four. 


The  analysis  occurs  in  several  distinct  steps,  each  corresponding  to  those  taken  in  the  exper¬ 
iments.  First,  a  temperature  history  is  applied  which  simulates  a  quench  of  the  composite  in  ice 
water.  In  computing  the  thermally-induced  residual  stresses,  the  material  is  taken  to  be  stress-free 
at  a  uniform  temperature  of  350° C.  Prior  studies  [7]  have  shown  negligible  dependence  of  computed 
residual  stresses  on  soiutionizing  temperatures  at  or  above  350° C  because  of  the  low  flow  strengths  of 
aluminum  alloys  at  these  temperatures.  Once  the  material  is  quenched,  the  mechanical  response  of 
the  composite  subject  to  unisixieJ  tension  or  compression  along  the  fiber  axis  is  simulated.  The  unit 
cell  is  then  unloaded  so  that  the  net  stresses  vanish  on  the  cell  walls,  and  is  subsequently  cooled  to 
liquid  nitrogen  temperature  (  — 196°C).  Finally,  the  composite  is  returned  to  room  temperature. 


Completing  the  equations  for  either  method  requires  specifying  the  appropriate  constitutive 
response  for  both  phases,  although  only  the  elastic  constants  are  required  for  the  approach  beised  on 
Eshelby  solutions  [3].  The  silicon  carbide  whiskers  are  taken  to  be  elastic  and  the  aluminum  matrix 
is  modelled  as  an  isotropically  hardening  elastic-viscoplastic  solid.  The  uniaxial  plastic  response  of 
the  matrix  is  given  by 


(2) 


with  a  reference  strain  rate  of  to  =  0.002s  \  a  reference  strain  of  to  =  0.004,  and  a  strain  rate 
hardening  exponent  of  m  =  0.004.  Poisson's  ratio  was  taken  to  be  i/  =  0.3,  and  the  coefficient  of 


thermal  expemsion  was  a  =  24  x  10“®  °K~^.  The  reference  strength  ao,  strain  hardening  exponent 
N,  and  Young’s  modulus  E  were  considered  to  be  functions  of  temperature.  Although  these  veJues 
were  known  at  lower  temperatures  from  ovu:  own  experiments,  material  data  was  not  available  for 
2009  aluminum  at  higher  temperatures.  The  properties  at  high  temperatures  were  assumed  to  be 
those  of  6061-T6  eduminum  [8],  which  has  roughily  the  same  flow  behavior  as  2009  aluminum  at  room 
temperature.  The  constitutive  relation  was  generalized  to  multi-axial  states  by  the  usual  procedure 
of  replacing  the  stress  and  plastic  strain  appearing  in  (2)  with  their  respective  von  Mises  equivcdents. 
The  property  values  of  silicon  carbide  were  taken  to  be  £  =  470  GPa,  u  =  0.19,  and  a  =  4.7  x  10“® 
[9]. 

Results  ztnd  Discussion 

The  stress  vs.  strain  curves  in  tension  and  in  compression  predicted  by  the  finite  element  simu¬ 
lations  of  a  unit  cell  reinforced  with  cylindrical  fibers  are  compared  with  experimented  data  in  Figure 
1.  The  test  specimens  exhibit  significantly  higher  flow  strengths  in  compression,  in  agreement  with 
other  published  reports  [10,11].  The  calculations  predict  nearly  identical  composite  flow  strengths 
by  3%  strain,  suggesting  that  damage  mechanisms  not  included  in  the  computations  may  cause  the 
observed  differences  between  the  tensile  and  compressive  stress-strain  behavior. 


Abs(True  Strain) 


Fig.  1  Comparisons  of  the  experimental  and  predicted  stress  vs.  strain  response  in  tension  and  in 
compression  of  a  2009  Al/  15  vol%  SiC  composite. 


The  predicted  levels  of  residual  axial  strain  in  the  matrix  as  a  function  of  composite  plastic 
deformation  for  the  three  different  models  considered  here  are  shown  in  Figure  2.  The  analytical 
method  of  Withers  [3]  et  al.  (here  assuming  no  thermal  residual  stresses)  gives  a  linear  dependence  of 
residual  strain  on  plastic  deformation.  For  comparative  purposes,  we  also  derived  solutions  for  a  simple 
radially-loaded  composite  spherical  shell  model,  assuming  an  incompressible,  elastic-perfectly  plastic 
matrix  reinforced  with  a  spherical  elastic  inclusion.  This  model  also  predicts  a  linear  relationship 
between  residual  strain  and  composite  plastic  strain.  In  contrast,  the  finite  element  solutions  predict 
a  non-linear  dependence.  Finite  element  ceilculations  with  ellipsoidal  inclusions  predict  roughly  50 
percent  lower  residual  elastic  strains  compared  to  those  obtained  with  a  cylindrical  fiber. 

The  residual  strains  measured  by  neutron  diffraction  are  compared  with  the  finite  element  pre¬ 
dictions  assuming  cylindriccil  reinforcements  in  Figure  3.  The  solutions  have  excellent  qualitative 
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Fig.  2  Comparisons  of  three  different  model  predictions  of  average  matrix  elastic  residual  strain  in 
the  axial  direction  vs.  imposed  composite  plastic  deformation.  Thermal  residual  stresses  were  not 
included  when  implementing  the  method  of  Withers  et  al.  [3]. 


agreement  with  the  experimental  curves,  and  very  good  quantitative  agreement  is  found  for  the 
residual  elastic  strains  in  the  aluminum  matrix.  Note  that  the  thermal  residual  strains  are  quickly 
dominated  by  the  residual  strains  generated  by  even  relatively  small  levels  of  composite  plastic  de¬ 
formation.  This  behavior  is  in  agreement  with  calculations  performed  in  [7],  which  demonstrated 
that  composites  with  and  without  thermally  induced  residual  stresses  included  in  the  aneJysis  had 
neeirly  identical  stress  states  after  small  amounts  of  deformation.  Also  note  that  the  residueil  elastic 
strains  axe  lower  in  the  fractured  specimen  (the  specimen  with  the  greatest  plastic  strain),  which 
would  suggest  that  the  composite  experienced  significant  damage  in  the  latter  stages  of  the  tensile 
test. 


Only  lattice  spacings  for  one  plane  of  silicon  carbide  could  be  measured  experimentally  in  each 
direction  because  of  the  alignment  of  whiskers  in  the  composite.  Although  the  trend  of  residual  axial 
strain  in  the  silicon  carbide  reinforcements  was  correctly  predicted  by  the  finite  element  calculations, 
the  magnitudes  of  residual  strain  were  consistently  underestimated.  In  contrast,  predicted  magnitudes 
of  residual  strain  in  the  transverse  direction  were  greater  them  the  measured  V2jues.  The  reasons  for 
these  discrepancies  has  not  been  conclusively  identified.  One  possibility  is  that  the  measurements 
of  whisker  axial  strains  involved  only  perfectly  aligned  whiskers,  and  therefore  may  not  have  been 
representative  of  the  average  state  of  the  reinforcements.  Another  possible  reason  could  be  attributed 
to  uncertainties  in  the  elastic  properties  of  the  whiskers.  In  addition  to  effects  of  elastic  anisotropy, 
reported  vsdues  of  elastic  moduli  for  ^  silicon  carbide  in  the  <  111  >  direction  range  from  352  GPa  [12] 
to  over  500  GPa  [13].  In  contrast,  the  aluminum  matrix  is  nearly  isotropic  and  its  elastic  constants 
are  well  characterized.  Finally,  the  silicon  carbide  diffraction  peaks  were  often  irregular  in  shape 
and  intensity,  and  consequently  lattice  spacings  obtained  using  the  curve-fitting  algorithms  exhibited 
greater  uncertainty  than  those  obtained  from  the  aluminum  diffraction  peaks.  These  irregularities 
may  be  associated  with  the  high  density  of  grown-in  defects  present  in  the  whiskers  [14]. 

Conclusions 


Mechanically  induced  residual  stresses  in  Al/SiC  composites  were  studied  with  analytical  eind 
numericsil  models,  and  with  neutron  diffraction  experiments  that  measured  average  residual  eleistic 
strains  in  the  fibers  and  matrix.  The  method  of  Withers  et  al.  [3]  predicted  a  linear  dependence 
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Fig.  3  Comparisons  of  average  residual  elastic  strains  vs.  composite  plastic  strain  as  measured  by 
neutron  diffraction  and  predicted  by  finite  element  calculations  with  cylindrical  reinforcements  for  fa^ 
the  zJuminum  matrix  in  the  axial  direction,  ^b)  the  aluminum  matrix  in  the  transverse  direction,  (c) 
the  silicon  carbide  whiskers  in  the  axial  direction,  and  (d)  the  silicon  czirbide  whiskers  in  the  transverse 
direction.  Included  in  the  aluminum  data  are  strains  measured  for  individual  lattice  planes  and  by 
Rietveld  profile  refinement  [6]. 


of  residual  elastic  strains  on  the  imposed  plastic  deformation,  in  disagreement  with  experimental 
results.  In  that  study,  the  authors  speculated  that  the  observed  non-linear  dependence  on  composite 
plastic  deformation  was  caused  by  room  temperature  diffusional  relaxation  of  the  aluminum  matrix. 
The  finite  element  calculations  suggest  that  this  is  unlikely,  since  good  quantitative  agreement  with 
experiments  was  observed  for  the  aluminum  matrix  without  any  diffusional  terms  included  in  the 
model.  Furthermore,  the  finite  element  solutions  also  suggest  that  inclusion  geometry  significantly 
affects  the  values  of  even  average  field  quantities. 

Only  qualitative  agreement  was  obtained  for  the  average  residual  strains  in  the  silicon  carbide 
whiskers.  Possible  reasons  for  the  discrepancies  include  fiber  misorientation,  poorly  cheuracterized 
elastic  moduli  for  the  silicon  carbide,  or  errors  induced  by  inconsistent  diffraction  peaks.  Because  the 
whiskers  are  essentially  rigid  in  comparison  to  the  plastic^ly  deforming  aluminum  alloy,  model  predic¬ 
tions  of  matrix  residu^  strains  are  relatively  insensitive  to  reinforcement  elastic  properties.  The  close 
correlation  of  the  numerical  model  with  neutron  diffraction  measurements  of  the  aluminum  matrix 
suggest  that  it  may  be  possible  to  predict  microstructural-level  stresses  with  unit  cell,  continuum- 
based  models  using  the  finite  element  method.  Future  work  will  address  residual  stresses  in  Al/SiC 
pcirticulate  composites  and  the  effects  of  more  complex  loading  paths. 
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Abstract 

The  precipitation  characteristics,  the  mecha¬ 
nisms  of  accelerated  aging,  and  the  variation 
of  uniaxial  tensile  stress-strain  behavior  in 
response  to  controlled  variations  in  matrix  micro- 
structure  were  investigated  for  a  2124  Al-SiC 
whisker  composite.  The  yield  strength  of  the  com¬ 
posite  was  found  to  he  independent  of  matrix  aging 
condition.  However,  the  overall  ductility  decreased 
monotonically  with  an  increase  in  aging  time.  A 
finite  element  analysis  of  the  constitutive  response 
of  the  composite  is  presented.  The  results  of  ihe.se 
calculations,  as  well  as  the  predictions  of  several 
models  for  composite  strengthening  available  in 
the  literature,  were  compared  with  the  experi¬ 
mental  results.  The  presence  ol  brittle  whiskers  m 
aluminum  leads  to  a  signilieani  build-up  of  hydro¬ 
static  stresses  in  the  matrix  during  plastu  detorma- 
non.  I  did  formation  in  the  matrix  of  the  composite 
as  well  as  at  the  wlusker-matrix  interlai  e  appears 
to  play  an  important  role  in  controlling  the  overall 
failure  mechanisms.  I  ransmissmn  eleitrof, 
microscopy  observations  ot  void  tonnatimi 
whisker  ends  are  described  tor  Kimpusin  \ihs  : 
mens  strained  in  tension  at  mom  tempemtioe  w;,.- 
at  300  X'.  A  detailed  discussion  ol  matrix  deimm..- 
tion  and  inierfacial  debonding  is  presented  in  a- 
attempt  to  identif\'  the  origins  oj  low  dm  tilits  m 
discantinuoiisly  reinforced  metal -ceramic  n'm 
poshes. 

I.  Introduction 

With  recent  advances  in  procewine  tech- 
nt'loev.  the  abilitv  to  produce,  in  econoniicallv 
feasible  quantities,  a  wide  r.inec  o|  met.il-niatnv 
composites  IS  increasine  ili.iin.iiK.illv  Metals 
reinforced  with  brittle  particles  olteii  reieited 
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to  as  discontinuously  reinforced  metal -matrix 
composites,  constitute  a  large  portion  of 
these  advanced  materials.  Particle-reinforced 
metal-matrix  composites  are  both  machinable 
and  workable  using  conventional  processing 
techniques,  and  the  manufacturing  of  these  mate¬ 
rials  can  easily  be  scaled  up  to  productiem  quan¬ 
tities.  t.xperimental  results  of  stress-strain 
relations  for  these  materials  have  been  reponed 
previously  il-.'sl.  although  the  effects  of  syste¬ 
matic  variations  in  the  matrix  microstructure  on 
the  overall  composite  properties  have  not  been 
examined.  Recent  work  has  demonstrated  that 
aging-induced  precipitation  and  microstructural 
evolution  in  the  matrix  of  whisker-reinforced 
composite  materials  can  be  significantly  accel¬ 
erated  compared  with  those  of  the  unreinforced 
matrix  material  subjected  to  identical  heat  treat¬ 
ments  o-s  These  results  imiicate  that  conven¬ 
tional  heat  treatments  for  unremtorced  com¬ 
mercial  aluminum  alloys  should  not  be  directly 
,i('plicd  to  the  microstructural  design  of  com- 
(s.'siiv  materials  although  this  is  often  the  practice 
'  u  III  \  thuiough  examination  of  micro¬ 
s’ -.ivtutai  cvoluiinn  111  the  matrix  of  the  composite 
.!  I’s  ettcus  .<n  I'verall  mechanical  behavior  is 
.  ’.v.i'  !.•  a  moK  complete  understanding  of  the 
.  Tis’itutivc  lespoiise  and  to  the  evaluation  of 
.  eqs. 'siiv  strengthening  theories.  The  low  duc- 
ev  exhibited  by  discontinuously  reinforced 
metal- matrix  composites  is  the  primary  obstacle 
preventing  their  introduction  into  many  structural 
applications.  \et  the  influence  of  controlled 
matrix  microstructural  v  ariations  on  the  strain-to- 
failure  of  whisker-reinforced  metal-matrix  com¬ 
posites  remains  poorlv  understood  This  report 
will  locus  on  cxpeiimental  observiitions  and  finite 
element  artaivso  oi  the  effects  of  mtitrix  micro- 
'iiusnii.il  van, Hams  on  the  unncxial  stress-strain 
icsjsmse  ot  a  discontinuouslv  reinlorcesi 
met.il  m.niiv  soniposiic  in  an  attempt  to  better 
unilcisiansi  the  lailure  mechanisms 
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The  material  selected  for  this  study  was  a 
2 1 24  Al- 1 5  wt.%(  1 3  vol.%  )SiC  whisker  reinforced 
metal-matrix  composite  obtained  from  ARCO 
Chemicals,  Greer  SC,  produced  using  powder 
metallurgy  techniques.  The  2 1 24  Al  matrix  mate¬ 
rial  had  a  nominal  composition,  in  weight  percent, 
of  4%Cu-2%Mg-().4%Mn,  with  the  balance  Al. 
For  comparison  where  appropriate,  an  unrein¬ 
forced  alloy  of  2 1 24  Al,  that  was  identically  pro¬ 
cessed  from  the  same  powder  batch  and 
.subjected  to  identical  solutionizing  and  aging 
treatments,  was  also  studied.  The  whiskers  used 
in  this  study  were  the  highest  quality  available 
and  were  classified  as  F'-y.  The  average  dimen¬ 
sions  of  the  whiskers  before  processing  were  0.5 
fim  in  diameter  and  25  //m  in  length.  Flowever. 
due  to  the  severe  deformation  involved  in  pro¬ 
cessing  of  the  composite  material,  significant 
breakage  of  the  whiskers  occurred  and  the  aver¬ 
age  length  of  the  whiskers  in  the  as-extruded 
material  was  about  2.5  yum.  A  systematic  micro- 
structural  characterization  and  aging  study  of  the 
same  extruded  bars  of  the  composite  and  control 
alloy  was  performed  by  Christman  and  Suresh  :6 
and  will  be  briefly  reviewed  in  a  later  section  the 
reader  is  referred  to  the  original  paper  for  a  more 
complete  discussion  i.  Furthermore,  composite 
damage  mechanisms  operative  during  uniaxial 
tensile  tests  at  high  temperatures  will  be  dis¬ 
cussed  and  compared  with  damage  mechanisms 
at  room  temperature. 

A  finite  element  unit  cell  model  was  used  to 
predict  the  uniaxial  stress  strain  response  of  the 
composite  material  as  a  function  of  matrix  mic- 
rosiructural  variations.  The  model  approximates 
the  stress -strain  re.sp(7nse  of  a  single,  rigid  cylin¬ 
drical  inclusion  in  an  elastic  power-law  hardening 
viscoplastic  matrix  material.  The  matrix  and 
whisker  parameters  were  taken  directly  from 
experiments  when  possible  i,  thereby  minimizing 
the  number  of  adjustable  variables.  The  experi¬ 
mental  observations  were  compared  with  the 
numerical  predictions  and  the  implications  of  the 
results  are  discussed  in  the  context  of  the 
mechanisms  responsible  for  the  high  values  of 
yield  strength  and  the  low  values  of  ductility  dis¬ 
played  by  the  material. 

2.  Microstructural  evolution 

Recently.  Christman  and  Suresh  6’’  conducted 
a  study  of  the  microstructural  evolution  and  aging 
kinetics  in  the  2124  Al-SiC  whisker  composite 


and  control  alloy  used  in  this  studv.  utili/ing 
analytical  transmission  electron  mieroscopv  tech¬ 
niques.  matrix  microhardness  measutenieritv  .nui 
conductivity  measurements.  I  he  aim  ol  the  work 
reported  in  ref.  6  was  to  examine  the  eltect'  >  >n  the 
matrix  of  the  composite  produced  In  iIk  imio 
duction  of  the  whiskers  in  controlled  mkf 
structures  by  a  comparison  ot  remior..cd  .md 
unreinforced  materials. 

The  composite  and  control  allov  were  icoei  ,  ed 
in  the  as-extruded  condition.  solution-iie.ii<.d  .it 
504  °C  for  4  h,  water  quenched,  and  subsequent K 
aged  at  177°C  for  various  times.  The  micro- 
structural  condition  of  the  matrix  material  was 
determined  by  microhardness  measurements 
iFig.  1  ai  and  by  analytical  electron  microscopy. 
U'ith  the  addition  of  the  SiC  whiskers,  the  peak 
aging  time  of  the  matrix  of  the  eompcwni.-  is  tirasn- 
calK  reduced  from  12  h  lor  the  coniiol  .illov  to 
about  3  h  .Another  inierestme  teaiurc  ot  ihi' 


o  « j - 

f  •  P«intofce<3  2i24 

.  A  Unremforcefl  2l24 


Fiu  1  I  he  vitriiition  nt  a  miciohardnes,  and  h  prccipiune 
size  with  ajnni;  lime  atiei  rel  n  . 


51 


graph  is  that  there  is  no  appreciable  vertical  shift 
in  the  microhardness  v  alues,  that  is  the  composite 
matrix  is  not  significantly  harder  than  the  unrein¬ 
forced  material.  This  result  suggests  that  it  would 
be  possible  to  use  the  bulk  mechanical  properties 
of  the  control  alloy  for  the  matrix  properties  of 
the  composite  in  theoretical  calculations  of  the 
constitutive  response  of  the  composite,  with 
proper  adjustments  for  changes  in  the  aging 
curves.  This  inference  has  a  strong  bearing  on  the 
evaluation  of  composite  strengthening  theories 
which  will  be  discussed  in  a  later  section. 

The  transmission  electron  microscopy  TEM 
results  showed  that  the  age-hardening  precipi¬ 
tates.  S',  evolve  into  corrugated  sheets  after 
nucleating  on  dislocations  and  growing  in  the 
three  preferential  (K)l  directions.  I  his  type  of 
precipitation  sequence  is  lypictil  ol  \I-C'u-Mg 
allov  systems  and  is  well  documented  ;n  the  liter¬ 
ature  1  1-1 -V-  There  is  little  preferential  precipi¬ 
tation  or  precipitate-free  zones  at  or  near  the 
whisker-matnx  interface  I  he  rate  of  nucleation 
of  the  precipitates,  however,  was  different 
between  the  reinforced  and  unreinforeed  materi¬ 
als  The  reinforced  material  showed  evidence  of 
S'  nucleation  after  less  than  1  h  of  artificial  tiging 
at  l'’7‘C'.  whereas  the  unreinforeed  material 
required  up  to  4  h  of  aging  at  this  temperature 
before  anv  S' could  be  detected.  Fig.  1  b  . 

I  he  accelerated  aging  of  the  composite  matrix 
material  is  aided  b\  a  decrease  in  the  incubation 
time  for  the  nucleation  of  age-hardening  precipi¬ 
tates.  I  he  large  thermal  contraction  mismatch 
between  the  aluminum  matrix  and  the  SiC 
whiskers  gives  rise  to  residual  stresses  upon  cool¬ 
ing  from  the  sulutionizing  temperature  Disloc¬ 
ations  are  punched  out  from  whisker  ends  during 
cooling  to  relieve  pan  of  the  residUiil  stresses  at 
the  interface.  For  2  1  24  .Al-  1  .^.2vol  '..SiC  compo¬ 
site.  the  dislocation  punching  distance  caleulaied 
from  ref.  14  is  sufficientlv  large  to  cover  the 
majority  approximatelv  iy'„  of  the  matrix  with 
excess  dislocations.  Transmission  electron  mic¬ 
roscopy  studies  by  Christman  and  Suresh 
indicate  that  the  dislocations  in  the  matrix  of  the 
composite  serve  as  nucleation  sites  for  the  fvtrma- 
tion  of  S  precipitates  during  the  aging  of  the 
composite.  This  facilitates  the  attainment  of  peak 
matrix  hardness  at  much  shorter  times  than  in  the 
control  alloy. 

The  theory  for  accelerated  aging  listed  above 
is  supported  by  a  wide  variety  of  evidence 
(Christman  and  Suresh  Ibj;; 


1 1 )  TEM  studies  of  the  composite  material  and 
unreinforeed  alloy  reveal  increased  dislocation 
density  in  the  matrix  of  the  composite. 

i2;  TEM  studies  also  reveal  that  the  strength¬ 
ening  precipitates  nucleate  preferentially  along 
dislocation  lines. 

i?'  Little  preferential  precipitation  or  precipi¬ 
tate  free  zone  at  or  near  the  whisker  interfaces 
was  observed. 

(4  Quantitative  analy  sis  of  the  growth  kinetics 
of  the  strengthening  precipitates  for  both  the  rein¬ 
forced  and  unreinforeed  materials  show  s  a  shift 
in  nucleation  time  but  not  in  growth  kinetics. 

'  Measurements  of  changes  in  electrical  con¬ 
ductivity  are  consistent  with  the  proposed  model. 

6  Cold  w  orking  of  the  unreinforeed  material 
also  produces  an  accelerated  aging  effect. 

7  Recent  in  siiu  TEM  experiments  have 
demonstrated  dislocation  punching  at  whisker 
ends  during  cooling  of  the  metal-matrix  com¬ 
posites  1  .s 

s  Calculations  ef  disk'cation  punching  dis¬ 
tance  agree  well  with  the  observed  dislocation 
distribution. 

Experimental  method 

1  he  composite  material  selected  for  this  inves¬ 
tigation  was  the  same  material  used  in  the  previ¬ 
ously  mentioned  characterization  study  b  and 
was  a  2124  powder  metallurgy  aluminum  alloy 
containing  1  ,s.2  vof’o  SiC  whiskers.  In  order  to 
form  a  basis  for  comparison,  a  control  alloy  with 
an  identical  processing  history  from  the  same 
powder  btitch  containing  no  reinforcement  was 
also  studied  The  material  was  received  as- 
extriulcvi.  in  bars  measuring  12.74  mmxl2'7 
mm  1  Son  mm  in  size.  T  he  heat  treatments  were 
designcil  to  achieve  identical  values  of  micro- 
hanlncss  m  the  unreinforeed  control  alloy  and 
in  the  matrix  of  the  reinforced  material.  The  heat 
treatment  schedules  were  chosen  to  obtain  two 
underaged  conditions,  the  peak  aged  condition, 
and  two  overaged  conditions  as  outlined  in  Tabic 
1.  .After  heat  treatment,  the  samples  were  stored 
in  a  commercial  freezer  in  order  to  prevent  any 
aging  at  room  temperature. 

Tensile  tests  were  performed  on  an  Instiam 
screw -driven  testing  machine  at  room  tempera¬ 
ture  Using  samples  conforming  to  .ASTNl  stan¬ 
dard  ES-S.s  for  sub-size  specimens.  The  samples 
were  machined  prior  to  heat  treatment  because 
machining  has  been  shown  to  alter  the  aging  state 


TABLE  1  All  samples  solutionized  504  C  for  4h  and  water 
quenched 
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of  the  material  [161.  The  strain  rate  applied  to  the 
samples  was  2.7x10  '  s  '  and  strain  was 
monitored  with  a  clip  gage  type  extensometer 
mounted  on  the  sample.  Load  and  strain  output 
were  recorded  directly  on  a  HP  model  7044 A 
XY  recorder  for  subsequent  analysis.  Vlicros- 
copic  examination  of  the  fracture  surfaces  was 
performed  with  an  AMR  model  lOOOA  scanning 
electron  microscope. 


4.  Finite  element  formulation  and  material 
model 

The  finite  element  analysis  is  based  on  a  con- 
vected  coordinate  Lagrangian  formulation  of  the 
field  equations  with  the  initial  unstressed  state 
taken  as  a  reference.  All  field  quantities  are  con¬ 
sidered  to  be  functions  of  convected  coordinates, 
y'.  which  serve  as  particle  labels,  and  time.  i.  This 
formulation  has  been  employed  extensively  in 
previous  finite  element  analyses,  see.  for 
examples,  refs.  17-19. 

/\ttention  is  confined  to  quasi-static  deforma¬ 
tions  and.  with  body  forces  neglected,  the  rate 
form  of  the  principal  of  \  irtual  work  is  w  ritten  as 


At  J  [f  dt  +  T  (i  ‘dn.,  ,  Id  I ' 

I 

=  AfJ  t  dudS-  J  T"dE„dV-\  T'dudS 


ill 


Here,  r  '  are  the  contravariant  components  of  the 
Kirchhoff  stress  x  =  J a.  with  a  being  the  Cauchv 
stress  on  the  deformed  convected  coordinate  net 
and  L  and  S  being  the  volume  and  surface, 
respectisely.  of  the  body  in  the  reference  con¬ 
figuration;  =  0  ,'0/  at  fixed  y  and  the  second 
term  on  the  right-hand  side  represents  an  equilib¬ 
rium  correction  term  that  is  used  in  the  numer¬ 
ical  procedure  to  reduce  drift  from  the 


Em.  1.  .4  schcnuuic  durnram  of  the  finiic  elcmeiil  mode!  t  he 
lenole  ihe  i  '  a.M'.  !■-  taken  to  be  ihe  eenierliiie  Hee.m^e 
of  '.ymmetr\  condilK'n's.  it  is  onK  neeessan  ii'  anaK/e  one 
quadrant  ol  the  unit  cell. 


equilibrium  path  due  to  the  discrete  time  step. 
The  nominal  traction  components.  /  .  and  the 
Lagrangian  strain  components.  H  ..  are  gi\  en  by 

=i  T'  +  T‘''n  i  v,  2 

=  !;  u, ,  +  II, ,  +  II , ''ll I,  ,  .7 

where  v  is  the  surface  normal  in  the  reference 
configuration,  ii,  are  the  components  of  the  dis¬ 
placement  vector  on  base  vectors  in  the  reference 
configuration,  and  i  ^  denotes  covariant  differen¬ 
tiation  in  the  reference  frame. 

A  cy  lindrical  coordinate  sy  stem  r.  6.z  <  is  used 
where  the  identifications  y'  =  r.y-  =  B  and  y  '  =  : 
are  made.  As  .shown  in  Fig.  2.  we  consider 
circular  cy  lindrical  fibers  of  radius  /y  embedded 
in  a  circular  cylindrical  cell  of  radius  R,.  and 
length  2L|,  and  with  an  initial  spacing  of  2l\. 
between  fiber  centers.  Attention  is  confined  to 
axisymmetric  deformations  so  that  all  field 
quantities  are  independent  of  B  and  within  each 
cell  symmetry  is  assumed  about  the  center  line. 
Furthermore,  the  circular  cylindrical  cell  sur¬ 
rounding  each  fiber  is  required  to  remain  a 
circular  cylinder  throughout  the  deformation  his¬ 
tory.  As  discussed  by  Tvergaard  |2(L  this 
axisymmetric  configuration  can  be  considered  an 


approximation  to  a  three-dimensional  array  of 
hexagonal  cylinders. 

The  boundary  conditions  for  the  axisymmetric 
region  analyzed  numerically  are: 

u''  =  0  /''  =  ()  onz  =  0  (4/ 

li  '  =  i’y  =  +  L'-j  t'  =  0  on  ;=/>,,  (5,i 

li '  =  t’l  7  '  =  0  on  r  =  /?,,  (6  i 

Here,  is  a  prescribed  constant  while  i\  is 
determined  from  the  condition  that  the  average 
lateral  traction  rate  vanishes,  i.e. 

h,. 

J  r'dc  =  0  on  r  =  7f,,.  i7) 


In  addition  to  the  boundary  conditions  !4i  to  {1 1, 
there  is  the  requirement  that  displacement  com¬ 
ponents  vanish  on  the  surface  of  the  rigid  fiber. 

In  some  calculations,  a  second  set  of  boundary 
conditions  was  employed  consisting  of  (4)  and 
iS).  but  with  if)}  and  (7)  replaced  by  7  '=()  on 
r  =  R„.  so  that  every  point  along  r  =  R,,  is  stress 
free.  Under  these  conditions  the  outer  sidewall  of 
the  cell  does  not  remain  straight  and  vertical.  Of 
course,  for  the  entire  tension  specimen.  7  '  =  0  on 
the  outer  boundary  of  the  specimen.  However, 
when,  as  here,  the  fibers  are  much  smaller  than 
the  specimen,  the  condition  that  the  circular 
cylindrical  cell  remain  a  circular  cylinder  comes 
from  enforcing  geometric  compatibility  (within 
the  axisymmetric  approximation!  for  a  uniform 
array  of  fibers  perfectly  aligned  with  the  tensile 
axis.  Relaxing  the  boundary  condition  i6j  and  ‘1  > 
permits,  in  a  highly  approximate  manner,  con¬ 
sequences  of  deviations  from  this  highly  con¬ 
strained  fiber  distribution  to  be  explored.  In 
subsequent  discussion,  predictions  based  on  '.6i 
and  ( 7 1  are  referred  to  as  results  "w  ith  constraint" 
and  predictions  based  on  7  '  =  ()  on  r=R„  are 
referred  to  as  results  "without  constraint". 

The  material  is  characterized  as  an  isotropi¬ 
cally  hardening  elastic  viscoplastic  solid  and  the 
total  rate  of  deformation.  D.  is  written  as  the  sum 
of  an  elastic  part.  D".  and  a  plastic  part.  D'’.  with 
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I  +  i'  . 
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(8) 


where  f  is  the  Jaumann  rate  of  Kirchhoff  stress.  I 
is  the  identity  tensor.  A ;  B  denotes  A  "77„.  e  is  the 


effective  .strain  rate.  E  is  the  Young's  modulus,  v 
is  Poisson's  ratio  and 


T'=T-3iT:l)l  d-  =  :T':T'  10 

/=  f„ld/gi  f ,ij‘  gi  f  I  =  o,,  f  /f„ -t-  1  T 

=  aJE  11 

Here.  f  =  /fd/  and  the  function  g  f;  represents 
the  effective  stress  I'.s.  effective  strain  response  in 
a  tensile  test  carried  out  at  a  strain  rate  such  that 
f  =  f,,.  Also.  o„  is  a  reference  strength  and  .V  and 
m  are  the  strain  hardening  exponent  and  strain 
rate  hardening  exponent,  respectively. 

Combining  eqns.  (8i  and  (0  and  inverting 
gives. 
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For  use  in  eqn.  i  1 standard  kinematic  rela¬ 
tions  are  used  to  express  eqn.  ( 1 2  as  a  relation 
between  the  contravariant  components  (on  the 
current  base  vectors ()  of  the  convected  rate  of 
Kirchhoff  stress,  f".  and  Lagrancian  strain  rate. 

The  deformation  history  is  calculated  in  a 
linear  incremental  manner  and.  in  order  to 
increase  the  stable  time  step,  the  rate  tangent 
modulus  method  of  Peirce  ei  al.  [2 1 J  is  used.  "This 
is  a  forward  gradient  method  based  on  an 
estimate  of  the  plastic  strain  rate  in  the  inter\al 
between  i  and  i+Ai.  The  incremental  boundary 
value  problem  is  solved  using  a  combined  finite 
element  Rayleigh-Ritz  method.  Tvergaard  [  1 9  . 


5.  Experimental  results 

Representatise  stress-strain  curves  for  the 
unreinforced  control  allo\  and  the  reinforced 
composite  material  are  shown  in  Fig.  3.  The  addi¬ 
tion  of  the  whiskers  causes  an  increase  in  the 
elastic  modulus.  ().2'’o  offset  yield  strength  and 
ultimate  strength,  and  a  decrease  in  the  strain  ti' 
failure  for  the  composite  material.  These  eftects 
are  well  documented  in  the  literature  [  1  M.  11'. 
What  is  more  surprising,  however,  are  the  effects 
of  aging  state  on  the  stress-strain  behavior  of  the 
composite  material.  Though  the  aging  curse  for 
the  unreinforced  control  alloy  isariation  of  0.2" o 
offset  yield  strength  as  a  function  of  aging  time 
exhibits  the  well  known  "bell-shaped"  curve,  the 
composite  yield  strength  is  independent  of  aging 
time  (Fig.  4ia  .  The  ultimate  strength  of  the 
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Pig.  3.  Reprosenuilivo  stress-strain  curves  for  the  composite 
matenal  lahelled  1 3"..  tmd  the  unreinlorced  control  allov 
labelled  ll".  .  The  increase  in  the  Youngs  modulus,  vield 
strength,  and  ultimate  strength  and  decrease  in  the  ductiliis 
for  the  composite  matenal  should  be  noted. 


Fig  4.  1  he  sanation  of  a  the  sield  strength  and  b  the 
ductiliis  ssith  aging  time 


compo.sile  (about  650  MPa.i  i.s  also  independent 
of  aging  time.  However,  the  strain-to-failure  of 
the  composite  is  a  strong  function  of  aging  time, 
decreasing  monotonically  from  about  5%  in  the 


as-quenched  condition  to  about  1.5‘’o  in  the 
severely  overaged  microstructure,  Fig.  4(bi. 

Microscopic  examination  of  the  fracture 
surfaces  revealed  a  dimpled  fracture  surface  for 
both  the  reinforced  and  unreinforced  materials. 
The  fracture  .surface  of  the  unreinforced  material 
contained  a  rather  even  distribution  of  large 
dimples  connected  by  sheets  of  smaller  dimples 
indicating  a  pattern  resulting  from  ductile  void 
growth.  coale.scence  and  failure.  Fig.  5.  The  frac¬ 
ture  surface  of  the  reinforced  material  contained 
only  .small  dimples  similar  to  the  small  dimples  on 
the  fracture  surface  of  the  unreinforced  material. 
Fig.  6(a).  The  aging  time  had  no  detectable  effect 
on  the  appearance  of  the  fracture  surface. 
Whisker  pullout  was  sometimes  observed; 
although,  the  pulled-out  whiskers  were  coated 
with  matrix  material.  Fig.  6(b  i.  This  demonstrates 
that  the  failure  of  the  composite  was  pre¬ 
dominantly  through  the  matrix  and  not  along  the 
matrix-reinforcement  interface.  Similar  conclu¬ 
sions  were  also  obtained  in  the  fracture  .studies  on 
SiC  whisker  and  particulate  reinforced  2XXX 
and  7XXX  alloys  |22]. 

6.  Observ  ations  of  void  formation  at  elevated 
temperatures 

One  mechanism  that  appears  to  contribute  to 
the  unusually  low  ductility  of  Al-SiC  whisker 
composites  involves  the  formation  of  voids  at 
fiber  ends.  Nutt  and  Duva  [23]  have  observed 
voids  forming  at  room  temperature  in  a  6061 
alloy  reinforced  with  SiC  whiskers  strained  to 
various  amounts.  Nutt  and  Needleman  124]  have 
developed  a  quantitative  description  of  the  void 
nucleation  process  at  fiber  ends  and  compared 
the  predicted  and  observed  modes  of  failure 
initiation.  While  further  studies  are  needed  to 
establish  the  link  between  void  formation  at 
whisker  ends  and  overall  failure  mechanisms  at 
room  temperature,  our  recent  work  indicates  that 
cavitation  at  whisker  ends  plavs  a  greater  role  in 
controlling  the  overall  constitutive  response  of 
Al-SiC  composites  at  elevated  temperatures  than 
at  room  temperature. 

The  initiation  of  voids  at  whisker  ends  was 
studied  in  the  same  21 24-SiC  whisker-reinforced 
composite  described  earlier.  When  the  test  tem¬ 
perature  is  increased,  the  tensile  ductilitv  of  the 
Al-SiC  composites  increases  significanth  ]25  . 
We  have  explored  this  temperature-dependent 
phenomenon  through  TE.M  observ  ations  of  spec- 


rii;.  5-  Fracture  >urtacc  morpholi>c\  ol  iIr  iimiitu. a , ; I  .  i.  .  •'!  the  i.'nipi  'itc  m.itc- 

C(mtn)l  allin  a  shnuiiic  the  Jivinbuiioii  .<1  ilti.  I.itci  itin.pi.  ,  ,,  j  ..  .  |  i_  '  be  iii'icil  \ 

and  b  vhowinj;  the  di'-lnbuluin  111  the 'titall  dtmplc'  ,  .a  •  ,■  ^  ,  ,■  ,  .'iix  ir.au  n.jl  i-  ^lii-ui:  in 

imens  tested  in  tensinn  at  3(Hi  C.  In  these  experi¬ 
ments.  \xe  ha\e  used  the  specimens  in  the  II  M  sp,.viinens  ucie  scelu'ned  trom  tensile 
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matrix  microstructure.  After  the  tensile  tests,  libei  eiiris  m  specimens  wnh  maeroscopic  plastic 


strains  ranging  from  3'’o  to  1  2. 5”,..  The  voids 
generally  formed  at  the  corner  on  the  fiber  end,  as 
shown  in  Fig.  Tiai.  and  intense  plastic  deforma¬ 
tion  occurred  in  the  vicinity  of  the  fiber  end.  as 
evidenced  by  the  high  dislocation  density  in  that 
region.  Glide  of  the  dislocations  was  inhibited  by 
the  coarse  age-hardening  precipitates  that  grew 


hit  ■  Void  niKk'.iimn  ;il  fiber  ends  m  2  I  24  .At -Si(  whisker 
eompi'siie  deliuiTied  ,il  .'Uli  r  in  ;i  tensile  sirnin  ol  11.12.  .i 
\i>id  ill  ei'rner  nl  liber  eiirl  surrounded  b\  dislueiiled  m;itri\. 
b  (  o.ileseenee  i>l  sokIs  toriiied  ,u  dittereril  siies  on  the  fiber 
end 


during  the  high  temperature  test,  although  the 
high  local  stresses  still  produced  large  plastic 
strains  in  this  region.  With  increasing  strain,  addi¬ 
tional  voids  formed  at  other  corner  sites  of  the 
same  fiber  end  Fig.  7  b".  and  eventually  the  voids 
coalesced  to  form  one  large  void  equal  m  width  to 
the  fiber  diameter. 

Alter  macroscopic  ten.sile  strains  of  12.."^",,, 
elongated  cavities  were  often  observed  at  fiber 
ends,  as  shown  in  Fig.  8  as  The  salient  features  of 
the  cavitation  mechanism  are:  1  i  the  cav  ity  width 
does  not  e.xceed  the  fiber  diameter.  2  there  is  no 
apparent  separation  along  the  fiber  side,  and  .s 
the  tip  of  the  cavitv  has  retained  the  same  basic 
shape  displaved  by  much  smaller  cavities 
associated  with  lower  strain  levels.  Furthermore, 
similar  elongated  cavities  were  observed  at  fiber 
breaks,  as  shcwvn  in  Fig.  8  b  .  These  observations 
indicate  that  large  amounts  of  shearing  at  or  near 
the  intetface  c:in  occur  after  decohesion  at  the 
fiber  end.  ana  that  the  load-bearing  capacity  <vf 
the  fibers  can  be  severelv  reduced. 

Patterns  of  void  evolution  observed  in  2  1  24 
Al-.SiC  composites  tested  at  .301)  "C  resemble  the 
observations  from  room-temperature  tests  in 
other  composite  systems  in  that  voids  nucleate 
near  the  corners  of  fiber  ends  and  subsequently 
coalesce  to  form  a  single  void.  However,  the  dis¬ 
tribution  of  voids  differs  in  that  vivids  nucleate 
throughout  the  gage  length  during  the  high-tem- 
perature  tests,  while  in  room-temperature  tests 
the  voids  are  confined  to  regions  verv  near  the 
fracture  surface.  Thus,  one  effect  of  temperature 
is  to  allow  the  damage  to  accumulate  in  the  com¬ 
posite  microstructure  prior  to  failure,  thereby 
increasing  the  ductilitv.  These  observations 
suggest  that  one  path  tcv  improving  the  room- 
temperature  ductilitv  of  Ai-SiC  composites 
would  be  to  identify  microstructural  parameters 
that  could  control  the  development  of  cav  itation. 
Toward  this  end.  calculations  based  on  an  inter¬ 
lace  decohesion  model  are  being  carried  out  to 
determine  the  effects  of  microstructural  param¬ 
eters  on  void  htrmation. 

7.  Discussion 

The  results  of  the  tensile  tests  conducted  in  this 
studv  are  useful  for  the  evnluaiion  of  composite 
strengthening  theories.  The  Modifievi  Shear  L  ag 
Morlel  .26  predicts  the  v  ield  strength  of  the  com¬ 
posite  materiiil  using  a  model  based  on  simple 
transfer  of  shear  across  the  vvhisker-mainv  inter- 


He  K  \oid  erowih  in  2H4  Al-Si(  whisker  Ciimposnc 
dclnrmcil  in  lensuin  ai  dmi  't'  a  I  hmeated  taxils  at  ihc 
fiber  end  siirrcunded  b\  dislDeallons  b  (Mindneal  earns  al 
the  fiber  break. 


fdce  vrith  an  additionfil  term  used  to  aeeount  Itn 
load  trtinsfer  across  the  ends  o(  the  w  hiskers.  'I  his 
theor\.  as  well  ;is  the  wideK  known  rule  ol 
mixtures  theory,  tire  not  in  eiiuilittitive  tiereement 
with  the  experimenttil  observations  ;is  they  would 
also  predict  th;it  the  yield  strength  ol  the  compo¬ 
site  mtiteruil  would  increase  with  an  increase  in 


the  yield  strength  of  the  matrix.  The  hnhanced 
Dislocation  Den'itv  Model  27]  attempts  to 
predict  the  yield  strength  of  discontinuously  rein¬ 
forced  metal-matrix  composites  based  on  a  mod¬ 
ified  shear  lag  model  taking  into  account  the 
differences  in  dislocation  densities  between  the 
reinforced  and  unreinforced  materials.  The  quali¬ 
tative  predictions  of  this  theory,  however,  are  in 
contradiction  with  the  experimental  observations 
of  the  present  study.  If  the  yield  strength  ol  the 
matrix  of  the  composite  is  significantiv  diflerent 
from  that  of  the  unreinforced  alloy,  then  one 
would  expect  to  see  this  difference  reflected  as  a 
uniform  increase  in  the  microhardness  values  for 
the  matrix  of  the  composite  material.  Figure  1  a 
indicates  that  the  matrix  of  the  composite  mate¬ 
rial  has  about  the  same  hardness  as  the  unrein¬ 
forced  alloy.  Also,  this  theory  wivuld  predict  that 
the  yield  strength  of  the  composite  material  will 
mirrvir  that  of  the  unreinforced  alloy  as  a  function 
of  aging  time,  although  this  is  not  observed 
experimentally.  Thus  the  Hnhanced  Dislocation 
Density  Model  is  not  in  qualitative  agreement 
with  the  experimental  results  and  therefore  is  of 
questionable  value  for  identifying  the  strengthen¬ 
ing  mechanisms  for  the  present  material. 

In  ref.  2S  a  self-consistent  formulation  for  the 
hardening  behav  ior  of  discontinuously  reinforced 
metal-matrix  composites  was  developed,  in  an 
attempt  to  predict  the  yield  strength  and  harden¬ 
ing  exponent.  The  analysis  in  ref.  2S  is  based  on  a 
path  independent  deformation  theorv  of  plasticilv 
anti  depends  on  the  use  of  llvushin's  theorem 
.2^:.  In  this  formulation  the  hardening  exponents 
of  the  reinforced  and  unreinforced  materials  are 
identical.  I'his  is  contrtidictorv  to  the  present 
experimental  observation  that  there  is  a  differ¬ 
ence  of  about  ;i  factor  of  three  between  the  har¬ 
dening  exponent  of  the  reinforced  and 
unreinforced  materials.  In  view  of  the  lack  of  an 
adequate,  analyticallv  tractable  model  for  the 
yield  and  htirdening  behavior  of  materials  rein¬ 
forced  with  discontinuous  fibers,  a  finite  element 
method  was  used  in  this  studv 

The  finite  element  method  possesses  scver.il 
inherent  advantages  over  manv  of  the  more  ana- 
Ivticiil  methods.  First,  geometrical  effects,  such  a-- 
corners  at  fiber  ends,  can  be  htindled  m  ;i  strtiight- 
forward  mtinner.  f  he  simpliiication  required  m 
manv  theoretical  devlopments  of  reducing  the 
inclusion  from  ii  evlinder  to  an  ellipsoiri  is  not 
necesstirv.  .Scconti.  the  calculation  of  stre^^  aiki 
strain  contour',  within  the  matrix  and  ;it  the 
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matrix-reinforcement  interface  are  possible. 
Third,  complete  stress-strain  responses,  as 
opposed  to  certain  portions  of  the  stress-strain 
curves,  are  predicted.  Therefore,  a  finite  element 
approach  is  appealing  for  this  application.  Here, 
an  a.\isymmetric  finite  element  unit  cell  model  is 
employed  where  the  aspect  ratio  of  the  whisker  is 
the  average  aspect  ratio  observed  experimentally 
and  the  aspect  ratio  of  the  unit  cell  is  taken  t(t  be 
approximately  equal  to  the  whisker  aspect  ratio. 
A  similar  finite  element  model,  but  using  a  full 
three-dimensional  cell  geometry,  is  described  in 
ref.  30. 

The  finite  element  method  predictions  for 
uniaxial  loading  are  qualitatively  in  agreement 
with  the  experimental  results.  Figure  doi;  shows 
the  predicted  stress-strain  curses  as  a  function  of 
the  N  olume  fraction  of  reinforcement  for  whisker- 
reinforced  2124  A1  in  an  underaged  state.  The 
effects  of  matrix  aging  condition  on  the 
stress-strain  response  are  plotted  in  Fig.  Oibi. 
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The  finite  element  model  correctly  predicts  that 
the  yield  strength  of  the  composite  material  is 
independent  of  matrix  variations  due  to  aging 
and  that  the  strain-hardening  exponent  of  the 
composite  be  quite  different  from  that  of  the 
unreinforced  material. 

The  quantitative  predictions  for  stress-strain 
lesponse.  however,  are  not  within  acceptable 
bounds  of  the  experimental  results.  Figure  lb 
shows  the  predicted  and  experimental 
stress-strain  curves  for  an  underaged  condition 
in  the  composite  material.  The  calculated  curves, 
with  a  traction-free  lateral  side  and  with  the  side- 
wall  constraint  representing  the  periodic  arras 
imposed,  bound  the  experimental  results, 
although,  they  are  not  sufficiently  close  to  be  of 
use. 

The  analysis  of  the  effects  of  the  sidewall  con¬ 
straint  on  the  resprmse  of  the  unit  cell  is  useful  in 
attempting  to  understand  the  underlying  strength¬ 
ening  mechanism  for  these  materials.  V\'ithout  the 
imposed  constraint,  there  is  little  increase  in  the 
yield  strength  of  the  composite  over  that  of  the 
unreinforced  control  alloy.  The  increase  in  yield 
strength  of  these  materials  appears  to  originate 
not  from  direct  load  transfer  to  the  reinforcement 
phase  but  from  the  constraint  of  geometrical 
compatibility. 

The  distribution  of  hydrostatic  stress  within 
the  unit  cell,  with  and  without  the  sidewall  con¬ 
straint.  elucidates  this  mechanism  further  hvr 
.dentical  geometries  and  aging  states  at  similar 
values  of  far-field  axial  strain.  Fig.  1  1 .  The  differ- 
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Hii;,  11.  Uisinhutum  ot  h\driistalic  in  GPa  wiihin  the 
unii  cell.  'Aiih  and  withinii  the  eide\^all  constraint.  1  he  values 
o(  tar  held  axial  strain  i'  almost  identical  in  both  cases.  The 
shaded  portion  ol  the  lijttire  represents  the  whisker. 


ence  in  far-field  axial  stre  is  about  300  MPa. 
which  IS  approximate!),  equal  to  the  difference  in 
the  \alues  for  hydrostatic  stress  within  the  upper 
region  of  the  cell.  As  far-field  strain  accumulates, 
the  sidewall  in  the  upper  portion  of  the  cell 
attempts  to  contract  at  a  greater  rate  than  in  the 
lower  portion  of  the  cell.  The  sidewall  constraint 
forces  the  lateral  side  to  remain  straight,  thereby 
generating  large  compressive  hydrostatic  stresses 
within  the  lower  portion  of  the  cell  and  large 
tensile  hydrostatic  stresses  within  the  upper 
portiem  of  the  cell.  It  is  these  tensile  hydrostatic- 
stresses  which  appear  to  increase  the  yield 
strength  of  the  composite  material  rather  than 
direct  load  transfer  to  the  reinforcement  phase. 

The  cell  model,  "with  constraint",  i.e.  based  on 
the  boundary  conditions  6-  and  7,.  represents 
an  ideal  regular  array  of  perfectly  aligned  fibers 
and  appears  to  overpredict  the  levels  of  hydro¬ 
static  stress  generated  within  the  material.  In 
certain  regions  of  the  cell,  the  tensile  hydrostatic 
stress  component  reaches  three  times  the  yield 
strength  of  the  matrix  material.  At  the  present 
time,  it  is  not  clear  that  such  high  levels  of  hydro¬ 
static  tension  are  actually  achieved  and  there  are 
a  variety  of  mechanisms  that  could  act  to  reduce 
the  hydrostatic  tension.  Tor  example,  deviations 


from  the  aligned  end-to-end  fiber  geometry 
assumed  in  the  calculations  would  be  expected  to 
lead  to  reductions  in  the  hydrostatic  tension 
levels  achieved.  Indeed,  preliminary  analysis  of 
this  effect  by  the  authors  indicate  significant 
reductions  are  so  attained. 

In  the  calculations  discussed  so  far  the  fiber 
has  been  taken  to  be  rigid.  A  few  calculations 
were  carried  out  accounting  for  fiber  elasticity.  In 
these  calculations  Young's  modulus  for  the  SiC 
fibers  was  taken  to  be  450  GPa  and  Poisson  s 
ratio.  v  =  0.1 7  [3 1 1.  Compared  with  the  rigid  fiber 
analysis,  there  is  a  reduction  in  the  Young's 
modulus  of  the  composite  of  8%- 10"...  Also,  the 
peak  hydrostatic  tension  is  reduced  by  5'’o-l()''.. 
when  fiber  elasticity  is  accounted  for.  These 
calculations  show  the  quantitative  magnitude  of 
effects  associated  with  fiber  elasticity  and  indicate 
that  in  the  Al-SiC  system  investigated  there  iv  no 
qualitative  difference  in  behavior  ansing  from  the 
assumption  of  rigid  fibers. 

The  present  finite  element  model  makes  no 
attempt  to  predict  ductility.  However,  the  predic¬ 
tions.  in  conjunction  with  the  experimental 
results,  provide  some  insight  into  possible  failure 
modes  at  room  temperature.  The  failure  of  the 
composite  appears  to  be  strongly  dependent  on 
the  matrix  material.  Recent  work  by  Christman 
and  Suresh  [32]  on  the  response  of  the  composite 
material  and  control  alloy  to  fatigue  loading 
demonstrated  that  under  these  conditions  the 
behavior  of  the  composite  is  controlled  by  the 
failure  of  the  matrix  material.  The  lack  of  exposed 
whiskers  on  the  fracture  surface  and  the  strong 
dependence  of  the  ductility  on  aging  raise  ques¬ 
tions  as  to  the  relative  roles  of  voids  at  the 
whisker-matrix  interface  and  those  within  the 
matrix  of  the  composite. 

The  large  amounts  of  local  deformation 
observed  on  the  fracture  surl'ace  in  the  form  of 
dimples  suggests  that  a  void  formation  mech¬ 
anism  is  present.  This  is  consistent  with  the  high 
levels  of  hydrostatic  stress  predicted  by  the  finite 
element  model.  Even  though  the  existence  of  the 
high  levels  of  hy  drostatic  tension  predicted  by  the 
finite  element  model  is  uncertain,  the  existence  of 
significant  tensile  hydrostatic  stresses  within  the 
matrix  material  appears  plausible.  The  generation 
of  these  localized  tensile  hydrostatic  stresses,  in 
the  matrix  slightly  above  the  whisker,  may  be  a 
driving  force  for  void  formation  in  the  matrix 
material  at  relatively  low  values  of  far-field  axial 
.strain. 


It  Nhould  be  stated  that  other  metal-matrix 
composite  materials  in  other  studies  have  dis¬ 
played  strikingly  difterent  tractographic  features 
and  failure  modes  (33.  34|.  The  dependence 
of  ductilit\  on  reinforcement  chemistry  and 
morphology,  matrix  chemistry  and  processing 
techniques  is  evident,  although  not  well  docu¬ 
mented  or  characterized.  The  need  for  systematic 
studies  of  the  dependence  of  ductility  in  metal- 
matrix  composite  systems  on  these  parameters 
is  clear.  A  basis  for  such  a  study  should  be  a 
thorough  in\estigation  of  ductility  in  a  single 
metal-matrix  composite  sy  stem. 


8.  Concluding  remarks 

In  this  study,  the  uniaxial  tensile  stress-strain 
response  of  a  2 1 24  Al-SiC  whisker  metal-matrix 
composite  was  investigated  as  a  function  of  con¬ 
trolled  matrix  microstructural  variations.  The 
(t.2"o  offset  yield  strength  of  the  composite 
material  was  independent  of  aging  time  despite  a 
significant  change  in  the  yield  strength  of  the 
unreinlorced  control  alloy.  Several  models  avail¬ 
able  in  the  literature  for  composite  strengthening 
were  reviewed  and  discussed  in  relation  to  the 
experimental  results  of  this  study  and  an  earlier 
microstructural  characterization  study  4>:.  .A 
finite  element  unit  cell  model  was  shown  to 
predict  the  correct  qualitative  trends,  although 
good  quantitative  accuracy  was  lacking.  The 
results  of  the  finite  element  model  vvere  discussed 
in  an  attempt  to  develop  a  better  understand¬ 
ing  of  the  strengthening  mechanisms  in  these 
materials. 

The  mechanisms  involved  in  the  tensile  quasi- 
static  failure  of  the  metal-matrix  composite  were 
investigated  at  room  temperature  and  at  30ii  (  . 
\'oid  formation  mechanisms  apireared  to  play  a 
dominant  role  in  the  failure  behavior  at  both  tem¬ 
peratures.  .Although  the  precise  role  of  the  voids 
in  the  fracture  process  remains  unclear  at  this 
time,  voids  nucleating  at  whisker  ends  appear  to 
be  more  predominant  at  elevated  temperatures 
than  at  room  temperature.  Further  investigation 
is  required  to  determine  the  dominant  charac¬ 
teristics  of  voids  formed  at  the  whisker-matrix 
interface  and  those  formed  within  the  matrix  of 
the  composite.  In  addition,  the  mechanisms  mllu- 
encing  the  strong  dependence  of  ductility  on 
aging  condition  remtiin  to  be  demonstrated. 
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Abstract 

The  effect  of  thermally  induced  residual  stresses  on  the  mechanical  propenies  and  ductility  of  Al-Si( 
composites  was  investigated  numerically.  The  predicted  behavior  in  unia.sial  loading  was  calculated  with 
and  without  the  inclusion  of  the  residual  stresses  which  result  from  the  mismatch  in  thermal  expansion 
between  aluminum  and  SiC.  In  this  analysis,  void  nucleation  by  interfacial  debonding  at  the  whiskers' 
ends  was  assumed  to  be  the  limiting  failure  mechanism.  Two  cases,  both  with  a  fiber  volume  fraction  of 
and  fiber  aspect  ratio  of  4.  but  with  different  fiber  spacings.  were  considered.  The  residual  stresses 
had  a  small  effect  on  the  predicted  ductility  of  the  composite,  even  when  a  relatively  weak  interface 
strength  was  assumed.  The  residual  stresses  are  shown  to  redistribute  as  interfacial  failure  is 
approached.  A  close  end-to-end  fiber  spacing  gives  a  greater  flow  strength  in  compression  than  in 
tension  and  the  residual  stresses  which  arise  during  thermomechanical  processing  tend  to  enhance  this 
effect. 


1.  Introduction 

Aluminum  alloys  reinforced  with  silicon 
carbide  whiskers  display  significant  increases  in 
stiffness,  tensile  strength,  and  creep  resistance 
when  compared  with  unreinforced  alloys  il-."^:. 
Unfortunately,  the  composites  also  exhibit  poor 
ductility,  which  has  thus  far  limited  their  applica¬ 
tion  as  a  structural  material.  A  better  understand¬ 
ing  of  the  underlying  mechanisms  affecting 
composite  properties  is  therefore  essential  if  the 
material  is  to  be  improved. 

Transmission  electron  microscopy  TENU  ob¬ 
servations  by  Nutt  and  Duva  |4]  revealed  that 
void  nucleation  at  the  w  hisker  ends  occurred  near 
the  .surface  of  fractured  specimens.  The  decohe¬ 
sion  of  the  matrix  from  the  whisker  end  at  low 
macroscopic  strains  appears  to  contribute  to  the 
poor  ductility  observed  in  Al-SiC  composites. 
Other  damage  mechanisms  observed  experi¬ 
mentally  include  void  nucleation  by  fiber  break¬ 
age  [5.  6]  and  void  growth  and  coalescence  in  the 
matrix  it.self  [.'i.  7], 

Nutt  and  Needleman  |8j  analyzed  Al-SiC  com¬ 
posites  within  the  framework  of  interface 
debonding  developed  by  Needleman  [9],  Further 
calculations  [lOj  predicted  a  transition  in  the 


mode  of  debonding  depending  on  fiber  spacing; 
debonding  initiated  at  the  whisker  comer  for 
close  side-to-side  spacings  of  fibers  but  began  in 
the  fiber  center  when  spacings  were  close  end-to- 
end.  Examination  of  a  fractured  Al-SiC  speci¬ 
men  16.  8]  using  TEM  tended  to  confirm  the 
results,  as  debonding  in  the  whisker  center  was 
generally  only  observed  when  fiber  spacing  was 
close  end-to-end.  The  strong  similarities  of 
observed  and  predicted  patterns  of  void  evolu¬ 
tion  suggested  that  it  may  be  possible  to  quantify 
interface  properties  by  comparing  experimental 
observations  and  predictions  based  on  finite 
element  calculations.  However,  the  calculations  in 
refs.  8  and  1(1  assumed  that  the  composite  was 
initially  stress  free. 

Arsenault  and  Taya  1 1 1 1  observed  yield 
strengths  higher  in  compression  than  in  tension 
for  a  606]  aluminum  alloy  reinforced  with  SiU 
whiskers.  In  addition,  they  modeled  the  matrix  as 
a  bilinear  solid  and  used  Eshelby's  solution  |12j 
for  an  ellipsoidal  inclusion  to  arrive  at  estimates 
of  residual  stresses  and  the  resulting  effect  on 
uniaxial  stress-strain  behavior.  Similar  experi¬ 
ments  and  analyses  have  also  been  performed  by 
Warner  and  Stobbs  [  1  V.  The  limitations  of  apply¬ 
ing  Eshelby's  solution  to  Al-SiC  composites 
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include  assuming  temperature  independent 
material  properties  and  ellipsoidal  whiskers. 
Approaching  the  problem  computationally 
allows  the  incorporation  of  more  realistic  matrix 
constitutive  behavior  and  fiber  geometries. 

Residual  stresses  are  inherent  in  Al-SiC  com¬ 
posites  owing  to  the  mismatch  in  thermal  expan¬ 
sion  between  aluminum  and  silicon  carbide.  The 
stresses  develop  upon  cooling  after  thermo¬ 
mechanical  processing  of  the  material.  In  pre¬ 
vious  work  [I4j,  we  inve.stigated  the  effects  of 
composite  microstructural  parameters,  such  as 
fiber  volume  fraction,  aspect  ratio,  and  spacing, 
on  the  resulting  residual  stresses.  In  many  cases, 
significant  compressive  stresses  were  calculated 
at  the  whisker  end.  Finite  element  residual  stress 
computations  have  also  been  performed  by 
Zivsicz  and  Parks  1151  for  continuous  fiber- 
reinforced  graphite-aluminum  composites.  The 
objectise  of  the  present  work  is  to  examine  the 
effects  of  the  thermally  induced  residual  stresses 
on  mechanical  properties  and  failure  of  Al-SiC 
whisker  composites  by  comparing  the  predicted 
response  when  residual  stresses  are  included  with 
the  response  when  the  material  is  initially  stress 
free. 

2.  Problem  formulation 

The  approach  used  here  is  to  model  the 
material  as  a  periodic  array  of  hexagonal  cylin¬ 
ders.  each  consisting  of  an  aluminum  alloy  with 
an  SiC  whisker  embedded  in  the  center.  The 
hexagonal  cylinder  can  then  be  approximated  as  a 
circular  cxiinder.  so  that  the  problem  is  then 
rendered  axissmmetric  see  Fig.  1  .  In  addition, 
because  of  the  symmetry,  onh  one  quadrant  in 
the  r-:  plane  need  be  analyzed. 

The  anal\  sis  can  be  thimghl  of  as  occurring  in 
two  distinct  steps.  First,  a  temperature  history 
simulating  a  quench  of  the  composite  is  applied  to 
the  unit  cell  depicted  in  Fig.  1.  The  cell  walls  are 
required  to  remain  straight  during  the  thermally 
induced  deformation,  which  maintains  compati- 
bilitx  with  surrounding  cells.  Additionally,  the 
average  stresses  along  the  walls  were  fixed  at 
zero,  simulating  the  condition  that  the  macro¬ 
scopic  stresses  during  quenching  are  zero.  Once 
the  residual  stress  state  in  the  composite  has  been 
determined,  the  mechanical  response  of  the  com¬ 
posite  sub|ect  to  uniaxial  tension  or  compression 
along  the  fiber  axis  is  simulated.  The  boundary 
conditions  on  the  cell  walls  are  that  thex  remain 


straight  during  loading  and  that  the  average 
stresses  on  the  side  cell  xvall  vanish.  I'he  bound¬ 
ary  conditions  in  both  cases  are  applied  by  use  of 
a  mixed  finite  element  and  Rayleigh-Ritz  method 
developed  by  Tvergaard  !16j.  The  mechanical 
response  of  the  composite  with  residual  stresses 
is  then  compared  with  the  response  of  the  mate¬ 
rial  that  is  initially  stress  free. 

The  fiber-matrix  interface  is  incorporated  into 
the  finite  element  equations  by  the  method  pre¬ 
sented  bx  Needleman  (yj.  We  note  that  although 
the  approach  is  purely  continuum,  the  finite 
element  formulation  of  the  field  equations  fulK 
accounts  for  matrix  plasticity,  thermal  expansion, 
general  temperature  dependence  of  material 
properties,  and  finite  geometry  changes. 

Completing  the  finite  element  equations 
requires  specifying  constitutixc  relations  for 
whisker,  matrix  and  interface.  The  SiC  whiskers 
are  taken  to  be  elastic  and  the  aluminum  matrix 
is  modeled  as  an  isotropically  hardening 
elastic-viscoplastic  solid.  The  uniaxial  plastic 
response  of  the  aluminum  matrix  is  gix  en  bx 
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The  properties  of  the  matrix  were  chosen  to  be 
representative  of  6061  aluminum  and  are  given 
by  a  reference  strain  rate  of  f„  =  0.002  s  a 
reference  .strain  of  f,,  =  0.004.  a  strain  rate  hard¬ 
ening  exponent  of  m  =  0.004.  Poisson  s  ratio  of 
V  =  0.3.  and  a  coefficient  of  thermal  expansion  of 
a  =  24xlo  '■  K  '.  The  reference  strength  a„. 
strain  hardening  exponent  .\'.  and  Young  s  modu¬ 
lus  £  were  considered  to  be  functions  of  tem¬ 
perature.  Using  material  data  at  various 
temperatures  1 1 7j.  the  properties  are  interpolated 
by  the  use  of  cubic  splines.  The  data  specified 
above  result  in  uniaxial  stress-strain  curves  at 
various  temperatures  shown  in  Fig.  2  in  ref.  14. 
The  constitutive  relation  is  generalized  to  multi- 
axial  states  by  the  usual  procedure  of  replacing 
the  stress  and  plastic  strain  appearing  in  eqn.  1 
with  their  respective  von  .Mises  equivalents.  .A 
complete  development  of  the  constitutive 
description  is  given  in  ref.  14.  The  property 
values  of  fi  SiC  are  given  bv  £  =  470  GPa. 
v  =  0.1 9.  and  a  =  4.7  x  lb  K  '  il^j. 

■A  potential  d  is  assumed  to  exist  for  the  inter¬ 
face.  which  is  taken  to  depend  only  on  the  jump 
in  displacements  across  the  fiber-matrix  inter¬ 
face.  From  the  potential,  the  interfacia'  iractions 
force  per  unit  initial  area  are  specified  by 


for  u„  <  d.  and  T„  =  /,  =  0  for  n,,  >  d.  where  n,,  and 
u.  represent  the  normal  and  tangential  difference 
in  displacements  across  the  interface  respectively. 
As  the  interface  separates,  the  magnitude  of  the 
nominal  tractions  increases,  reaches  a  maximum, 
and  finally  falls  to  zero  when  complete  decohe¬ 
sion  occurs.  From  eqn.  2  the  wxirk  of  separation 
at  !<,  =  ()  is  given  by 


The  interface  can  thus  he  characterized  by  the 


strength  o,,,,,.  the  work  of  separation  o  ,  .  and  the 
shear  parameter  y.  so  that  the  value  of  6  o  given 
through  eqn.  3  *.  The  parameter  6  is  a  pheno¬ 
menological  parameter  which  characterize'-  the 
ductility  of  the  separation,  and  has  dimensions  ot 
length. 

In  the  cases  studied  here,  the  interface  was 
characterized  by  =  3a,,.  where  a,,  is  the  room 
temperature  yield  strength  of  aluminum. 
<i  =  0.(l05  ,wm.  and  7=1.0.  Based  on  FF.M 
observations  of  voids  and  calculations  performed 
in  ref.  X.  the  value  of  the  interface  strength  used 
here  probably  represents  a  weak  Al-SiU  intei- 
face  Comparisons  of  the  most  commonly  ob¬ 
served  void  shapes  and  nucleation  sites  in  ret.  .s 
suggest  that  interface  strengths  ot  .sa,  -6a,.  are- 
more  likely  for  aluminum  bonded  to  SiC.  Low 
values  for  interface  strengths  were  used  in  the 
analy  sis  to  highlight  the  effect  of  residual  stresses 
and  because  material  damage  is  likely  to  initiate 
at  poorly  bonded  intertaees. 

3.  .Numerical  results 

In  computing  the  residual  stresses,  the  material 
was  taken  to  be  stress  free  at  a  uniform  tempera¬ 
ture  of  /,,.  z\  temperature  history  ot  the  form 
7  /  =  /,,  exp  -y.lf't'  with  t  in  seconds  vvas 
applied  uniformly  over  the  entire  cy  linder.  Such  a 
history  is  representative  (4  the  temperatures 
experienced  by  the  material  when  undergc'ing  a 
quench  in  ice  water.  The  spatial  uniformity  of 
temperature  in  the  unit  cell  is  justified  by  calcu¬ 
lations  in  ref.  14.  which  shovved  that  tor  realistic 
quenching  rates,  plastic  deformation  in  the  matrix 
IS  unlikely  to  induce  localized  heating.  Calcul¬ 
ations  were  performed  with  initial  temperatures 
of  either  /,,=  35'.'°C  or  /,,  =  .^6(i'C.  although 
material  properties  were  held  constant  bevond 
35(CC  due  to  lack  of  available  data.  The  depen¬ 
dence  of  the  resulting  residual  stresses  on  the 
assumed  temperature  of  the  initial  stress-free 
state  was  found  to  be  small.  This  is  anticiptitcd 
because  of  the  low  flow  strength  of  aluminum  at 
high  temperatures.  The  residual  stresses  used  in 
examining  the  mechanical  behavior  of  the  com 
posite  were  calculated  assuming  the  material  wav 
quenched  from  7,,  =  3.''('  T'. 

One  result  of  the  study  in  ref.  14  was  that  the 
residual  stresses  at  the  whisker  end  were  found  to 
depend  primanly  on  the  side-to-side  spacing  of 
the  fibers.  In  this  investigation,  two  different 
spacings  of  the  fibers,  henceforth  referred  to  as 
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cases  (a)  and  (b).  were  considered,  both  with  a 
fiber  aspect  ratio  of  4  and  a  fiber  volume  fraction 
of  20%.  Case  (  a  '  has  a  cell  aspect  ratio  of  4,  giving 
a  ratio  of  fiber  radius  to  cell  radius  of  r,,/ 
=  0.585,  while  case  (b)  was  given  a  cell  a.spect 
ratio  of  2,  yielding  a  ratio  of  fiber  radius  to  cell 
radius  of  r,,//?,,  =  0.463.  Decreasing  corre¬ 
sponds  to  increasing  the  side-to-side  spacing  of 
fibers  and  to  decreasing  the  end-to-end  spacing 
of  fibers.  The  calculated  average  compressive 
residual  stress  at  the  whisker  end  for  case  (a  i  was 
175  MPa  while  case  (bi  had  a  value  of  401  MPa. 
It  is  interesting  to  note  that  identical  calculations 
performed  with  a  rigid  interface  gave  average 
compressive  stresses  of  176  MPa  and  454  MPa 
for  cases  (ai  and  ibt  respectively.  Apparently,  the 
added  compliance  of  the  interlace  model  can 
affect  the  predicted  results  if  the  residual  stresses 
are  high  enough.  This  is  particularly  true  near  the 
fiber  corner  where  very  high  residual  stresses  are 
calculated  around  750  MPa  for  case  b  with  a 
rigid  interface  . 

.r  L  Intcrfacuil  diXi^licsion 

Figure  2  showk  the  unit  cell  stress  is.  strain 
response  of  cane  (.t  i.  the  case  with  the  closer  side- 
to-side  fiber  spacing,  with  and  without  the  initial 
residual  stresses  included  in  the  analysis.  One  can 
interpret  the  curves  as  the  predicted  macroscopic 
stress  r.v.  strain  response  of  a  composite  if 
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Fig  ;  Stress  is  strain  response  for  case  a  .  ihe  case  wiih  a 
fiber  aspect  ratio  of  4.  a  cell  aspeci  ratio  of  4  and  a  fiber 
volume  fraction  of  211"..  closer  side-lo-side  spacing  .  with 
and  without  residual  stresses  included 


debonding  from  every  whisker  end  occurred 
simultaneously.  The  maximum  stress  point  on  the 
curves  corresponds  approximateK  to  the  point  at 
which  void  nucleation  begins.  In  the  actual 
material,  one  could  imagine  damage  propagating 
through  the  composite  by  initial  debonding  at  a 
specific  whisker,  followed  by  load  transfer  to 
surrounding  regions  and  nucleation  of  new  v  oids. 
In  Fig.  2.  the  residual  stresses  have  little  efiect  on 
either  the  average  strain  to  debonding  or  the 
average  stress  at  which  debonding  initiates. 

Figure  3  displays  contours  of  axial  stresses 
near  the  maximum  load  point  for  case  a  w  ith  the 
residual  stresses  included  and  with  the  composite 
initially  stress  free.  The  two  stress  states  arc 
nearly  identical,  implying  that  the  thermally 
induced  residual  stresses  are  almost  completclv 
redistributed  by  the  time  the  composite  reaches 
maximum  load.  To  further  illustrate  this.  Fig.  4 
plots  the  average  normal  stress  at  the  fiber  end  i  v 
the  average  applied  stress.  Given  the  magnitude 
of  the  initial  residual  stresses,  the  predicted  insen¬ 
sitivity  of  the  overall  response  to  the  presence  of 
residual  stresses  is  somewhat  surpnsing.  Figure  4 
shows  that  the  difference  in  interface  normal 
stress  decreases  as  the  applied  load  increases,  and 
becomes  insignificant  near  the  point  of  maximum 
applied  .stress.  The  deformed  finite  element  mesh 
for  case  a'  at  an  average  strain  of  f  =  0.(i323  is 
shown  in  Fig.  5. 


W/  THERMAL  STRESSES  WVO  THERMAL  STRESSES 

Fig.  3.  Axial  strevses  for  case  a  .  the  case  with  closer  vrU- 
lo-side  spacing,  near  maximum  load  o,  =2i'.  soih 

initial  residual  stresses  and  with  the  composite  initialK  stress 
free 


Fii;,  -I  A\crai!i.'  '•trc'.sc^  at  tht  intcrtaci'  is.  ascraec  applied  Fii;.  Stress  n  strain  response  lor  ease  h  .  the  case  ssith  a 

stress  lor  case  a  .  the  case  ss  ith  closer  side-to-side  spacini;-  fiber  aspect  ratio  ol  4.  a  cell  aspect  ratio  ol  2  and  a  liber 

solume  Iraclion  ol  20"  ■  closer  end-lo-cnd  spacine  .  ssith  and 
ssithout  residual  stresses  included 


to-side  spacine.  at  an  aserace  strain  of  r  =11.1)32.''  shossine 

complete  debondine  at  the  sshisker  end.  Kic.  2  Aseratte  stresses  at  the  interdace  es.  aseraee  applied 

stress  for  case  b  ,  the  case  ssith  closer  end-to-end  spacing 

The  calculated  stress  vs.  strain  response  for 

case  !h!,  the  case  with  the  closer  end-to-end  fiber  Given  that  the  magnitude  of  residual  stresses  at 

spacing,  with  and  without  residual  stresses  the  interface  is  much  higher  in  case  (b)  than  in 

included  is  shown  in  Fig.  6.  Here  the  residual  case  lai,  a  greater  effect  was  expected.  Figure  7 

stresses  have  a  somewhat  greater  effect  on  the  displays  the  average  stresses  at  the  interface  i  v. 

predicted  mechanical  response  of  the  composite.  the  average  applied  stress.  Again,  as  w  ith  case  i  a  . 
Although  the  stress  at  which  debonding  occurs  is  as  the  average  stress  on  the  u  it  cell  approaches  a 

largely  unaffected  by  residual  stresses,  a  small  maximum,  the  difference  in  interface  stresses 

increase  in  the  strain  to  debonding  is  observed.  becomes  minimal. 
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Fig.  8  Stress  \s  strain  response  in  tension  and  compression 
for  case  a  .  the  case  with  a  cell  aspect  ratio  of  4  closer  side- 
to-side  spacing  :  a  with  residual  stresses  Included,  and  h 
with  the  composite  initialK  stress  free. 


S.2.  Mechanical  properties 

Figure.s  8  and  9  .show  the  effect  of  residual 
stresses  on  the  uniaxial  response  in  tension  and 
compression.  In  these  calculations,  the  focus  is  on 
the  mechanical  response  of  the  composite  and 
perfect  bonding  between  fiber  and  matrix  is 
assumed. 

Figure.s  8(a)  and  8(b)  respectively,  show  the 
uniaxial  response  in  tension  and  compression  of 


Fig.  SI.  Stress  is  strain  response  in  tension  and  compression 
for  case  h  .  the  case  with  a  cell  aspect  ratio  of  3  closer  end- 
to-end  spacing  :  a  with  residual  stresses  included,  and  b 
with  the  composite  iniiialK  stress  tree 


case  (a),  the  ca.se  with  the  closer  side-to-side  fiber 
spacing,  with  and  without  residual  stresses 
included  in  the  computations.  The  unit  cell  actu¬ 
ally  yields  at  a  lower  .stre.ss  in  compression  than  in 
tension  for  case  (a),  although  by  an  average  strain 
of  2%  the  two  curves  are  approaching  one 
another.  The  reduction  of  yield  strength  in  com¬ 
pression  is  not  particularly  surprising  since  sig¬ 
nificant  portions  of  the  matrix,  most  notably  near 
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the  fiber  comer,  have  compressive  residual  stress 
.states.  The  response  in  tension  and  compression 
is  essentially  identical  for  case  (a)  without  residual 
stresses  included. 

Figures  9(a)  and  9(b)  show  the  computed 
results  for  case  (b).  the  case  with  the  closer  end- 
to-end  fiber  spacing.  It  should  be  noted  that  here 
the  residual  stresses  have  the  oppo.site  effect, 
increasing  the  yield  strength  in  compression. 
Also,  the  yield  strengths  are  asymmetric  for  case 
(b),  even  when  no  residual  stresses  are  included. 
Figure  10  shows  the  large  plastic  strains  that 
develop  between  closely  spaced  fiber  ends.  This 
figure  illustrates  that,  at  a  given  level  of  average 
strain  in  compression,  there  is  greater  plas;;. 
deformation  in  the  matri.x  between  the  fiber  ends 
than  when  the  same  average  strain  level  is  applieo 
in  tension.  The  absolute  value  of  the  hydrostatic 
stress  between  the  fiber  ends  is  also  greater  in 
compression  than  in  tension.  For  close  end-to- 
end  fiber  spacings,  finite  geometry  change  effects, 
coupled  with  the  constraint  of  the  perfectly 
bonded  fibers  on  the  matrbc,  result  in  greater 
plastic  flow  between  the  whisfcer  ends  during 
compressive  loading. 


4.  Conclusions 

The  effect  of  thermally  induced  residual 
stresses  on  void  nucleation  in  Al-SiC  composites 
has  been  investigated  by  use  of  the  finite  ele¬ 
ment  method.  The  calculations  suggest  that  the 
residual  stresses  have  only  a  small  effect  on  either 
the  predicted  strength  or  ductility  of  the  material, 
at  least  when  ductility  is  limited  by  interfacial 
decohesion.  For  the  .-"satcrial  parameters  con¬ 
sidered  here,  the  residual  stresses  essentially 
become  fully  redistributed  as  maximum  load  is 
approached.  It  is  worth  bearing  in  mind  that  the 
cohesive  strength  used  in  the  calculations  prob¬ 
ably  represents  a  weak  Al-SiC  interface.  Some 
calculations  were  performed  with  a  stronger 
interface  and  the  residual  stresses  had  an  even 
smaller  effect  on  the  predicted  ductility  of  the 
composite  than  in  the  cases  shown  here.  This 
suggests  the  possibility  that  residual  stresses  may 
not  directly  have  a  significant  effect  on  ductility  of 
materials  where  the  predominant  damage 
mechanism  is  a  plastically  deforming  matrix  sepa¬ 
rating  from  hard  inclusions  ie.g.  spheroidized 
steels,  as  well  as  metal  matrix  composites), 
although  further  calculations  and  experiments 


Fig  in.  Effective  plastic  strain  contours  for  case  h  .  the  case 
with  a  cell  aspect  ratio  of  2  closer  end-to-end  spacing  .  and 
with  the  material  initialK  stress  free:  a  at  an  average  applied 
strain  =  0.(1273;  ib  at  an  average  applied  strain 

=  -0.0260.  Note  the  larger  straining  between  whisker 
ends  in  compression. 

would  be  required  for  confirmation.  However, 
residual  stres.ses  can  alter  the  microstructural 
development  of  matrix  alloys  [19j.  which  can  in 
turn  affect  the  ductility  of  the  composite.  Further¬ 
more,  the  residual  stresses  can  still  have  an 
important  effect  on  mechanical  properties,  as 
illustrated  in  Fig.  9. 

Another  important  damage  mechanism  ob¬ 
served  in  Al-SiC  composites  is  void  nucleation 
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by  fiber  breakage.  Figures  4  and  5  show  that  the 
initially  large  residual  compressive  axial  stresses 
in  the  fiber  ia  calculated  maximum  of  470  MPa 
and  980  MPa  for  cases  (ai  and  (b)  respectively) 
are  relieved  during  tensile  loading.  While  fiber 
breakage  is  not  incorporated  into  the  present 
model,  the  results  suggest  that  whether  or  not  the 
residual  stresses  significantly  inhibit  whisker 
failure  would  depend  on  whether  whisker  failure 
occurs  before  or  after  the  residua)  stresses  are 
relieved.  It  should  also  be  noted  that  there  may  be 
a  significant  effect  of  residual  stresses  on  defor¬ 
mation  and  failu’'e  mechanisms  not  considered 
here. 

Comparisons  of  the  mechanical  behavior  of 
the  composite  in  compression  and  tension  show 
that  the  effect  of  residual  stresses  is  sensitive  to 
fiber  spacing.  The  initial  flow  strength  in  com¬ 
pression  was  lower  than  in  tension  for  case  (ai. 
the  case  with  the  closer  side-to-side  fiber  spacing, 
although  perhaps  more  importantly,  the  overall 
flow  behavior  in  tension  and  compression  was 
quite  similar.  In  contrast,  case  (b  i.  the  case  with  the 
closer  end-to-end  fiber  spacing,  showed  signifi¬ 
cantly  higher  flow  strengths  in  compression. 
However,  it  should  be  recalled  that  even  with  'bi 
taken  to  be  initially  stress  free,  greater  flow 
strengths  were  predicted  in  compression  than  in 
tension.  In  fact.  Fig.  9  shows  that  at  an  applied 
strain  of  approximately  the  How  strength 
difference  between  tension  and  compression  is 
comparable,  with  or  without  residual  stresses 
included  in  the  analysis. 

1  he  present  results  indicate  that  close  end-to- 
end  fiber  spacings  promote  a  greater  flow 
strength  in  compression  than  in  tension.  The 
residual  stresses  which  arise  during  thermo¬ 
mechanical  processing  tend  to  enhance  this  effect. 
The  net  contribution  of  residual  stresses  to  the 
observed  difference  between  tensile  and  com¬ 
pressive  flow  strengths  depends  on  the  distribu¬ 
tion  of  fiber  spacings. 
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CREEP  DAMAGE  MECHANISMS  IN  CERAMIC  COMPOSITES 
S.  R.  Nutt  and  P.  Lipetzky 

Division  of  Engineering,  Bro\wn  University,  Providence,  Rl  02912 

Whisker-reinforced  ceramics  show  improvements  in  a  variety  of  mechanical  properties,  and  because  of  the 
inherent  thermal  and  chemical  stability  of  ceramics,  it  is  natural  to  consider  issues  related  to  mechanical  behavior 
at  high  temperatures.  Implementation  of  ceramic  materials  in  high-temperature  structural  applications  is  often 
limited  by  creep,  and  recent  studies  on  AljOg-SiC^  have  shown  that  creep  resistance  of  ceramics  can  be 

enhanced  by  whisker  reinforcement.  These  studies  demonstrated  that  creep  response  is  strongly  affected  by  the 
presence  of  the  reinforcing  phase  and  the  properties  of  the  whisker-matrix  interface.  In  addition,  oxidation 
reactions  tiiat  occurred  during  creep  degraded  interface  properties,  facilitated  microstructural  damage,  and 
accelerated  creep.  The  objective  of  the  present  study  was  to  determine  the  effects  of  test  ambients  on  creep 
damage  mechanisms  in  AlgOg-SiC^  composites  and  to  correlate  these  observations  with  compressive  creep 

response  under  different  stress-temperature  Conditions. 

Hot-pressed  SiC  whisker-reinforced  alumina  (from  Greenleaf  Corporation)  was  selected  for  the  present  study. 
The  composite  material  was  characterized  by  an  average  grain  size  of  1-2  microns  and  a  whisker  volume  fraction 
of  1/3.  The  composites  were  almost  completely  free  of  glass  phase,  at  least  within  the  detectability  limits  of  weak 
beam  dark  field  and  phase  contrast  imaging  (Fig.  1).  Material  was  machined  into  cylindrical  samples  for 
compressive  creep  experiments  which  were  conducted  in  air  and  in  nitrogen  at  1200-1400’C  under  applied 
stresses  of  25-250  MPa.  Creep  rates  were  monitored  until  a  nominal  constant  strain  rate  was  reached  and 
maintained.  Applied  loads  were  incremented  for  fixed  temperatures,  although  the  total  accumulated  strain  was 
kept  below  3%.  Creep  experiments  were  interrupted  during  secondary  creep  and  cooled  under  load  to  retain 
microstructures.  Thin  foil  specimens  were  sectioned  from  the  center  of  the  samples  and  prepared  by  ion  milling. 

Data  from  creep  experiments  performed  in  nitrogen  and  in  air  are  summarized  in  Fig.  2.  The  effect  of  the 
oxidizing  ambient  was  to  increase  the  creep  rate  by  factors  of  3-5  for  temperatures  of  1300-1400’C,  respectively. 
In  addition,  the  stress  dependence  of  the  creep  rate  increased  as  the  temperature  increased  above  1200'C.  The 
most  prevalent  damage  mechanism  observed  in  all  crept  samples  was  void  nucleation  at  grain  boundary-interface 
(GBI)  junctions,  an  example  of  which  is  shown  in  Figure  3.  Decohesion  occurred  after  accumulation  of 
aluminosilicate  glass  at  GBI  junctions  during  creep,  resulting  in  glass-lined  cavities.  In  the  non-oxidizing  ambient, 
the  aluminosilicate  phase  originated  from  the  small  amounts  of  residual  glass  films  present  in  the  as-fabricated 
material  (Fig,  1),  which  redistributed  during  creep.  In  the  presence  of  air,  the  amount  of  glass  phase  appearing  in 
the  composites  was  increased,  a  phenomenon  attributed  to  thermal  oxidation  of  SiC  near  the  specimen  surface. 
Silica  glass  produced  by  oxidation  during  creep  penetrated  the  composite  along  interfaces  and  grain  boundaries, 
accumulating  at  GBI  junctions  and  other  stress  concentrations.  An  example  oi  the  initial  stages  of  thermal 
oxidation  is  shown  in  Figure  4,  a  phase  contrast  image  of  a  near-surface  SiC  whisker  after  prolonged  creep  at 
1400'C  in  air.  Graphitic  cart)on  and  silica  glass  have  formed  at  the  interface,  although  the  carbon  subsequently 
disappears  via  formation  of  CO.  In  all  of  the  crept  samples  observed,  the  dislocation  density  was  negligible, 
leading  to  the  conclusion  that  the  composites  had  deformed  by  diffusional  creep.  However,  at  creep 
temperatures  above  1200‘C,  sliding  of  grain  boundaries  and  interfaces  became  increasingly  unaccommodated, 
especially  in  the  oxidizing  ambient.  Glassy  films  arising  from  thermal  oxidation  of  SiC  weakened  internal 
boundaries  in  the  material,  causing  extensive  microstructural  damage  (cavitation)  and  accelerated  creep.* 

‘Work  supported  by  the  Office  of  Naval  Research  through  Contract  #N00014-86-K-0125. 
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FIG.  1. -Residual  glass  film,  approximately  1  nm  thick,  at  whisker  interface  in  as-fabricated  composite. 

FIG.  2. -Steady-state  strain  rates  vs.  applied  stress  for  creep  experiments  in  nitrogen  and  in  air  at  1200-1400‘C. 
FIG.  3.-Glass-lined  cavity  at  GBI  junction  resulting  from  creep  in  air  at  1400’C. 

FIG.  4.-Thermal  oxidation  of  near-surface  SiC  whisker  during  creep  at  1400‘C  in  air,  resulting  in  silica  glass  fG) 
and  graphitic  carbon  (C). 
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